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ABSTRACT 

We report the detection of 48 millisecond pulsars (MSPs) out of 75 observed thus far using the LOw-Frequency ARray (LOFAR) 
in the frequency range 110-188 MHz. We have also detected three MSPs out of nine observed in the frequency range 38-77 MHz. 
This is the largest sample of MSPs ever observed at these low frequencies, and half of the detected MSPs were observed for the first 
time at frequencies below 200 MHz. We present the average pulse profiles of the detected MSPs, their effective pulse widths, and 
flux densities and compare these with higher observing frequencies. The flux-calibrated, multifrequency LOFAR pulse profiles are 
publicly available via the European Pulsar Network Database of Pulsar Profiles. We also present average values of dispersion measures 
(DM) and discuss DM and profile variations. About 35% of the MSPs show strong narrow profiles, another 25% exhibit scattered 
profiles, and the rest are only weakly detected. A qualitative comparison of the LOFAR MSP profiles with those at higher radio 
frequencies shows constant separation between profile components. Similarly, the profile widths are consistent with those observed 
at higher frequencies, unless scattering dominates at the lowest frequencies. This is very different from what is observed for normal 
pulsars and suggests a compact emission region in the MSP magnetosphere. The amplitude ratio of the profile components, on the 
other hand, can dramatically change towards low frequencies, often with the trailing component becoming dominant. As previously 
demonstrated this can be caused by aberration and retardation. This data set enables high-precision studies of pulse profile evolution 
with frequency, dispersion, Faraday rotation, and scattering in the interstellar medium. Characterising and correcting these systematic 
effects may improve pulsar-timing precision at higher observing frequencies, where pulsar timing array projects aim to directly detect 
gravitational waves. 
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1. Introduction 


Radio millisecond pulsars (MSPs) are commonly held to have 
been spun up by mass transfer from a binary companion (Alpar 
St al. 1982; Radhakrishnan &. Srinivasan 1981) and are evolu- 
tionarily linked with X-ray binary systems (Bhattacharya &. van 
den Heuvel 1991; Archibald et al. 2009; Papitto et al. 2013; Stap- 
pers et al. 2014). Their observed properties are in many respects 
quite different from those of the long-period, ’normal’ pulsars. 
Both the spin periods, P, and period derivatives, P, of MSPs 
[P ~ 1-lOms, P ~ 10^'^ss^') are a few orders of magnitude 
lower than for normal pulsars (P ~ 1 s, P ~ 10 s s '), and 
thus the inferred values of surface magnetic fields for MSPs are 
three to four orders of magnitude lower than the typical value of 
10'^ G for normal pulsars. Despite these differences, MSP stud¬ 
ies at 1.4, 2.7, and 4.9GHz by Kramer et al. (1998, 1999) and 
Xilouris et al. (1998) suggest that the radio emission mecha¬ 
nisms for normal pulsars and MSPs are essentially the same - 
even though their magnetospheres differ in size by two to three 
orders of magnitude. 


PSR J0437-4715 by Jenet et al. (1998), who find that its indi¬ 
vidual pulses are very similar to those of slowly rotating pulsars. 
Finally, Edwards & Stappers (2003) came to the same conclu¬ 
sion in their study of individual pulses from several pulsars with 
the Westerbork Synthesis Radio Telescope (WSRT) at several 
frequencies between 328 and 2240 MHz. 

Even though the radio emission mechanism of both MSPs 
and normal pulsars seems to be the same, there are a number of 
differences in the observed emission properties. Multifrequency 
studies above 400 MHz generally' show remarkably little fre¬ 
quency evolution in the average profiles of MSPs; in particular, 
the width of components and their separation typically remain al¬ 
most constant (Kramer et al. 1998, 1999). This is different from 
what is generally observed in normal pulsars (see e.g. Thorsett 
1991; Xilouris et al. 1996; Pilia et al. 2015), and may result 
from a small radio emission region in the very compact mag¬ 
netospheres of MSPs. Kramer et al. (1998) also show that MSPs 
are somewhat less luminous and less efficient radio emitters than 
normal pulsars, with isolated MSPs being even less luminous 
than their binary counterparts. These findings are in accordance 


This conclusion is mainly based on the similarity of their 
spectra and profile complexity. The same result was also ob¬ 
tained by Toscano et al. (1998) for a sample of MSPs in the 
southern hemisphere at frequencies between 0.4 and 1.6 GHz. 
This is also confirmed by single-pulse studies at 1.4 GHz of 


' There are notable exceptions, such as PSRs J2145-0750 and 
J0751-1-1807, that show large frequency evolution in the amplitude ratio 
of their profile components (Kramer et al. 1999; Kuzmin & Losovsky 
1996), and a few other MSPs have new profile components appearing at 
higher/lower frequencies. 
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with results from recent works by Levin et al. (2013) and Bur- 
gay et al. (2013). All these results suggest that MSP magneto¬ 
spheres might not simply represent scaled versions of the mag¬ 
netospheres of normal pulsars (Xilouris et al. 1998), and the dif¬ 
ferences could be caused by the different evolutionary histories 
of normal pulsars and MSPs. 

MSP studies at low radio frequencies, i.e. below < 300 MHz, 
are valuable, as they expand our ability to compare flux-density 
spectra, polarisation and profile evolution to those of normal 
pulsars and probe the compactness of MSP emission regions 
(Cordes 1978a). Normal pulsars often show a spectral turnover 
between 100 and 250MHz (e.g. Sieber 1973; Kuzmin et al. 
1978; Izvekova et al. 1981; Malofeev 2000; Malofeev et al. 
2000), raising the question of whether MSPs behave similarly. A 
study of 30 MSPs with narrow bandwidths at frequencies of 102 
and 111 MHz by Kuzmin & Losovsky (2001) concluded that the 
spectra of the MSPs they studied mostly do not show any low- 
frequency turnover. Only in the case of PSR J1012-I-5307 did 
they find a possible turnover near 100 MHz. Kuzmin & Losovsky 
(1999) analysed the frequency dependence of profile widths of 
12 MSPs between 102 and 1400 MHz. They found that profile 
widths remained nearly constant within this frequency range, 
and in the case of PSR J2145-0750, the component separa¬ 
tion became even smaller at 102 MHz than at higher frequen¬ 
cies (Kuzmin 8c Losovsky 1996). Both results provide further 
evidence for the compactness of the emission region of MSPs. 

Although Kuzmin & Losovsky (2001, and references 
therein) provided a first insight into low-frequency MSP pro¬ 
files and flux densities, their data suffered from poor time res¬ 
olution (only 0.128-0.64ms). This was mainly due to the ne¬ 
cessity to form many narrow frequency channels to minimise 
dispersion smearing. Observations of MSPs at such low fre¬ 
quencies indeed present a challenge because of the large pulse 
broadening caused by both scattering and dispersion in the in¬ 
terstellar medium (ISM). Dispersion in the cold ionised plasma 
of the ISM scales as Though computationally intensive, it 
can be fully accounted for by coherent dedispersion (Hankins 
1971; Hankins 8c Rickett 1975), unlike scattering, which also 
has an even stronger frequency dependence (scales as for 
a Kolmogorov spectrum of electron density inhomogeneities in 
the ISM). Because the pulses of the MSPs have much shorter 
durations than those of normal pulsars, they can be completely 
scattered out, preventing detection. 

On average, pulsars have a steep negative spectral index of 
-1.4 (Bates et al. 2013), and MSP spectra are found to be con¬ 
sistent with the normal pulsar population (Kramer et al. 1998). 
However, the background sky temperature of the Galactic syn¬ 
chrotron radiation is also frequency dependent, and this depen¬ 
dence is steeper Lawson et al. 1987). Thus, even if one 

assumes no turnover in the spectra of MSPs above 100 MHz, the 
majority of the MSPs, especially those located along the Galac¬ 
tic plane, will be more difficult to detect below 250 MHz than at 
higher frequencies. 

In the past decade, higher-quality, low-frequency MSP ob¬ 
servations have been enabled by a new generation of pulsar 
backends that provide wide-bandwidth, coherent dedispersion 
(sometimes in real-time). For example, Stappers et al. (2008) 
used WSRT and the PuMall pulsar backend (Karuppusamy et al. 
2008) in the frequency range of 115-175 MHz to coherently 
dedisperse the data to a final time resolution of 25.6-51.2 qs. 
They detected eight out of 14 MSPs observed and did not And 
any clear relationship between the detectability of a source and 
any combination of its period, DM or flux density at higher fre¬ 
quencies. Nevertheless, the majority of the non-detections were 


pulsars with lower flux densities. Stappers et al. (2008) did not 
detect two MSPs, PSRs J1024-0719 and J1713-H0747, that were 
detected by Kuzmin & Losovsky (2001)^. They did not And 
any obvious reason for these non-detections. The authors sug¬ 
gested scattering to be the main reason for all of their six non¬ 
detections. For PSR J2145-0750 Stappers et al. (2008) also 
found the reduced separation between the peaks of the compo¬ 
nents with regard to the high frequencies and even found this 
separation to be smaller than in Kuzmin & Losovsky (1996), 
which can be attributed to residual dispersive smearing. How¬ 
ever, Stappers et al. (2008) cast doubt on whether these compo¬ 
nents are the same as at high frequencies, because the leading 
component at 150 MHz more resembles the trailing component 
in the high-frequency profile. Recently, Dowell et al. (2013) also 
found the dominant component at low frequencies to be narrower 
than reported by Kuzmin & Losovsky (1996) in their observa¬ 
tions with the Long Wavelength Array (LWA) at 37-81 MHz. 

Since the work of Stappers et al. (2008) and until the re¬ 
cent study of PSR J2145-0750 by Dowell et al. (2013), there 
have been no other MSP studies at frequencies below 200 MHz 
that have provided high-quality profiles at these low frequencies. 
At the same time, the number of known MSPs in the Galac¬ 
tic field (many at distances < 2kpc, and at high Galactic lat¬ 
itudes) has substantially grown, owing to many detections in 
Fenni/LAT unassociated sources (see e.g. Keith et al. 2011; Ran¬ 
som et al. 2011; Kerr et al. 2012; Barr et al. 2013b). About 75 
MSPs were discovered in the Galactic field between the discov¬ 
ery of the first MSP (PSR B1937 -f 21; Backer et al. 1982) and 
the Fermi launch in 2008. In the last five years we have expe¬ 
rienced an MSP renaissance as this number has almost tripled 
to 218 (ATNF catalogue^; Manchester et al. 2005, and ongo¬ 
ing pulsar surveys). This is mainly a result of Fermi discoveries, 
but also includes a number of discoveries by ongoing pulsar sur¬ 
veys. These include the Green Bank Telescope (GBT) 350-MHz 
Drift-scan (Boyles et al. 2013; Lynch et al. 2013) and Green- 
Bank Northern Celestial Cap (GBNCC; Stovall et al. 2014a) 
surveys, the Pulsar Arecibo L-band Feed Array (PALFA) sur¬ 
vey (Nice et al. 2013; Crawford et al. 2012; Deneva et al. 2012), 
and the Parkes/Effelsberg High Time Resolution Universe Sur¬ 
vey (HTRU; e.g. Keith et al. 2010, 2012; Barr et al. 2013a). 
Therefore, there is also a significantly larger sample of MSPs 
that can be observed for the first time below 200 MHz. 

The LOw-Frequency ARray (LOFAR) is a digital aperture 
array radio telescope operating from 10-240MHz - i.e. the 
lowest four octaves of the ‘radio window’ observable from the 
Earth’s surface (see van Haarlem et al. 2013 for a description of 
LOEAR). It is an interferometric array of dipole antenna stations 
distributed over Europe, with the Core of 24 stations located 
within an area approximately 2 km across near the Dutch vil¬ 
lage of Exloo. Compared to previous low-frequency radio tele¬ 
scopes, LOEAR offers many advantages in its observing capa¬ 
bilities (Stappers et al. 2011). Among them are the large frac¬ 
tional bandwidth (80 MHz of bandwidth in the 10-90 or 110- 
250 MHz range), sensitivity that is at least 20 times better than 
the low-frequency frontends at WSRT (Stappers et al. 2008), and 
an unrestricted ability to track sources from rise to set, thus accu¬ 
mulating a large number of pulses in a single observing session. 
The latter is very valuable as many of the old low-frequency ar- 


^ Pulsar J0218+4232 was also not detected by Stappers et al. (2008), 
though Kuzmin & Losovsky (2001) reported it to be detected and pro¬ 
vided its flux density estimates, but did not present its profile. 

^ http://WWW.atnf.csiro.au/people/pulsar/psrcat/ 
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Table 1. Standard setup for the LOFAR observations of MSPs. 


Array 

Core 

stations 

Bits 

per 

sample 

Frequency 

range 

(MHz) 

Sub-bands" 

Sub-band" 

width 

(kHz) 

Sampling 

time 

(M-s) 

Observing 

time 

(min) 

Data 

products 

HBA Dual 

23^^ 

8 

110-188 

51-450 

195.312 

5.12 

20 

XXYY" 


Notes. The second polyphase filter (van Haarlem et al. 2013) was skipped in the online processing, thus the number of sub-bands and channels 
was the same, and we will refer lo the original sub-bands as channels in the rest of the paper. The core station CS013 has an incorrect orientation 
of dipoles with respect to the rest of the array and was always excluded in our observalions. Typically Iwo-lhree other core stations were not 
available for differenl reasons, and we used only 20-21 HBA stations out of 24 total. On a few occasions we only had 18-19 HBA slalions 
available. Raw complex vollages of Iwo linear polarisalions. 


rays, such as the Large Phased Array in Pushchino (Russia), are 
transit instruments. 

In this paper we present the results of the first exploratory ob¬ 
servations of MSPs with LOFAR. In Sect. 2 we describe the LO¬ 
FAR observing setup and data analysis. In Sect. 3 we present the 
results of our MSP observations, give lists of detected and non- 
detected pulsars, show folded pulse profiles and provide mea¬ 
surements of flux densities at 110-188 MHz. In Sect. 4 we dis¬ 
cuss the detectability of MSPs, provide new measurements of 
their DMs and compare these with previously published values. 
We also discuss DM and profile variations of MSPs and qualita¬ 
tively compare profiles for several pulsars with profiles at higher 
frequencies. We give a summary in Sect. 5. 


2. Observations and data reduction 

The observations presented in this paper were taken between 19 
December 2012 and 3 November 2014 with the LOFAR Core 
stations using the High-Band Antennas (HBAs) in the frequency 
range 110-188MHz and the Low-Band Antennas (LBAs) from 
10-88MHz (see van Haarlem et al. 2013; Stappers et al. 2011 
for more information on the LOFAR observing system and set¬ 
up). A single clock system allows all Core stations to be com¬ 
bined into one or multiple coherent tied-array beams, and an 
eight-bit sampling mode allows 80 MHz of instantaneous band¬ 
width (limited by filters that avoid the FM frequency range from 
90-110MHz). This provides a factor ~ 5 increase in sensitivity 
compared with early LOFAR pulsar observations that were re¬ 
stricted to using the six stations on the ‘Superterp’ and 48 MHz 
of bandwidth (e.g. Kondratiev et al. 2013). 

2.1. Data acquisition 

After coherent addition of the data streams from all avail¬ 
able Core stations (typically 20-23 stations), the raw complex- 
voltage (CV) data were recorded at 5.12 ps. The 78-MHz band¬ 
width is split into 400 channels of 195.312kHz each. The typical 
length of our observations was 20 min, although subsequent ob¬ 
servations for monitoring and timing purposes were adjusted de¬ 
pending on the pulsar’s signal-to-noise ratio (S/N) in the LOFAR 
band. 

The standard observing setup is summarised in Table 1. Dur¬ 
ing the observations the data from the stations were streamed for 
further online processing (beamforming, etc.) to the Blue Gene/P 
(BG/P) supercomputer"^ in Groningen. The output data products 


^ After April 18, 2014 (MJD 56765) the new LOFAR GPU 
correlator/beam-former Cobalt was used for MSP observations. 


(in our case the raw CV data) were written in HDF5^ format to 
the CEP2 cluster® for further offline processing. 

2.2. Puisar pipeiine 

After each observation, we ran the LOFAR PULsar Pipeline 
(PULP), a Python-based suite of scripts that provides basic of¬ 
fline pulsar processing (dedispersion and folding) for standard 
observing modes including CV data. PULP collects informa¬ 
tion about the observation itself and location of the data on the 
CEP2 nodes and then executes data reduction tools from the 
PRESTO^ (Ransom 2001), dspsr* * (van Straten & Bailes 2011) 
and PSRCHIVE® (Hotan et al. 2004) software suites, and creates 
a number of summary diagnostic plots and PSRFITS*® (Hotan 
et al. 2004) data files. 

To process CV data we used dspsr to read raw LOEAR data 
via the LOEAR Direct Access Library (DAL) interface. Each 
channel was coherently dedispersed and folded with a pulsar 
ephemeris. Where possible, we used the most recent published 
timing models for folding, but in some cases we used updated 
ephemerides from ongoing timing campaigns that are not yet 
published. We set the length of the folding sub-integrations to 
be 5 s and the number of pulse profile bins equal to the closest 
and lowest power of two for P/At (where P is the MSP period, 
and At - 5.12 ps is the sampling interval), but not higher than 
1024 bins per period. The data were combined in frequency us¬ 
ing psradd. In early observations, we used paz -r to perform 
zapping of a radio frequency interference (REI) with a median 
smoothed automatic channel zapping algorithm. In later observa¬ 
tions, we used the clean.py tool from CoastGuard" (Lazarus 
et al., in prep). This REI excision routine with the -F surgical 
option provided much better excision, and for this publication 
all pipeline products were reprocessed using this tool. The data 
were then dedispersed, i.e. appropriate time delays between fre¬ 
quency channels (note that each channel was already coherently 
dedispersed with dspsr) were applied, using pam. In the last 
step we ran the pdmp tool that enables the determination of a 
best DM and P by optimising the S/N of the profile. The output 
data files were further updated with these best DM and P values. 
Diagnostic plots were then produced to enable quick inspection 
of the data quality. 


^ Hierarchical Data Format 5, http: //www.hdfgroup. org/HDF5/ 

* LOFAR’s 2nd CEntral Processing (CEP2) computer cluster. It will 
be replaced with the next generation CEP4 cluster by the start of 2016. 
’ http: //WWW .CV.nrao.edu/~sransom/presto/ 

* http://dspsr.sourceforge.net 

^ http://psrchive.sourceforge.net 
http: //WWW .atnf.csiro.au/research/pulsar/index.html? 
n=Main.Psrfits 

* * https://github.com/plazar/coast_guard 
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3. Results 

At present we have observed 75 Galactic field MSPs with the 
LOFAR HBAs, and detected 48 of them. This is the largest sam¬ 
ple of MSPs ever observed/detected at these low frequencies. Of 
these sources, 25 are recent MSP discoveries and have never be¬ 
fore been detected at 110-188MHz. We discuss the detectabil¬ 
ity of MSPs in more detail in Sect. 4.1. The lists of detected 
and non-detected MSPs are given in Tables 2 and 3, respec¬ 
tively. For the detected MSPs we list the profile significance of 
the best individual observation and its S/N, observing epoch/ID, 
and integration length, Fobs- For both detected and non-detected 
MSPs we also list whether they were previously observed with 
the WSRT (Stappers et al. 2008; Archibald et al. 2009 for PSR 
J1023H-0038) or by the Large Phased Array (BSA) in Pushchino 
at 102/111 MHz (Kuzmin & Losovsky 2001). From the 27 MSPs 
that we did not detect, only one, PSR J0613-0200, was observed 
by Stappers et al. (2008) and they also did not detect it. Three 
out of these 27 have flux density measurements by Kuzmin & 
Losovsky (2001) (including PSR J0613-0200). However, they 
did not publish the profiles of these MSPs and their flux density 
estimates have very large uncertainties. 

3.1. Pulse profiles 

We present the best individual profiles for all 48 detected MSPs 
in Fig. 1. They are the sum of ~ 3 x 10^ spin periods, ranging 
from about 3 x 10^ for the 59.8-ms PSR J2235-H1506 to about 
1.3 X 10® periods for the 1-h observation of PSR J0337-I-1715 
(the pulsar triple system; Ransom et al. 2014). 

Many of the MSPs are detected with high S/N. For example, 
the black widow pulsar J1810-I-1744, discovered in the 350-MHz 
GBT Survey of faint Fermi y-ray sources (Hessels et al. 2011), 
was detected with S/N = 65 using a 20-min LOFAR observa¬ 
tion. Furthermore, the eclipsing pulsar J1816-I-4510, discovered 
in the GBNCC pulsar survey (Stovall et al. 2014a), was detected 
with S/N = 159 using a sum of four 5-min LOFAR observa¬ 
tions. LOFAR profiles for four PSRs J0214-H5222, J0636H-5129, 
J0645H-5158, and J1816-I-4510 were also presented in Stovall 
et al. (2014a) together with the GBT profiles at 350, 820, and 
1500 MHz. About 35% of the detected MSPs show strong nar¬ 
row profiles, while only 25% of MSPs have scattered profiles 
with a clear exponential tail. The profiles of the remaining MSPs 
(~ 40%) are weak, which could be due to the intrinsic spec¬ 
trum, and/or due to scattering, where the exponential scattering 
tail may approach or exceed a pulse period. Later in this Sec¬ 
tion we discuss the HBA detections of a few individual MSPs, 
and in Sect. 4.5 we discuss the profiles in general and make a 
comparison to high frequencies. Flux-calibrated, multifrequency 
LOFAR profiles from this paper are available via the European 
Pulsar Network (EPN) Database of Pulsar Profiles'^ (Lorimer 
et al. 1998). 

3.1.1. Folding and Worb/c caution 

Eor binary MSPs, the low observing frequencies and large frac¬ 
tional bandwidth of the LOEAR data mean that special care is 
needed when folding the data. This is because, in standard anal¬ 
yses, phase-folding is done based on a timing model evaluated at 
the central frequency of the band. If the dispersion delay across 
the band is large enough for the orbital phase to change signif¬ 
icantly during this time, then data can be assigned an incorrect 


http://www.epta.eu.org/epndb/ 


orbital phase. In such a case, the monochromatic timing model 
will not be able to correctly predict the pulse period across the 
entire bandwidth, but will have a frequency-dependent bias that 
changes with orbital phase (Ransom et al. 2004). Given that the 
error introduced in the pulse times-of-arrival (TOAs) is of the 
order forb/c x AfDM, where AfoM is the dispersion delay between 
the central frequency and the lowest frequency in the band, it 
is clear that this effect can most easily be countered by split¬ 
ting up the band into narrower sub-bands, and then applying 
the timing model separately to each of these. Alternatively, two- 
dimensional phase predictors can be employed. 

Without using two-dimensional phase predictors the actual 
error in TOAs will depend on the orbital phase of the pulsar and 
the projection of the orbit on the plane of the sky, and is largest 
for edge-on orbits when the pulsar is at periastron - i.e. when 
the pulsar’s radial acceleration is strongest. In general, the ef¬ 
fect is larger for higher-DM pulsars in compact (short orbital 
period) binary systems. In our MSP sample, the double pul¬ 
sar J0737-3039A has the largest value of i)orb/c ~ 10“^. Eor 
pulsars J0621-H1002, J2145-0750, J0751-bl807, J1023-b0038, 
J1022-H1001, J2215-H5135, and J1816-H4510, it is about (1- 
2)xl0 and for other binaries it is of the order of 10 ® or less. 
Using a single ‘polyco’ file (a file generated by the pulsar timing 
program Tempo'^ which contains a polynomial description of 
the apparent pulse phase as a function of time) applicable to the 
central frequency of the total band, would result in TOA errors 
of 8.4ms or 0.37 of the spin phase at the lowest observed fre¬ 
quency of 1 lOMHz for the pulsar J0737-3039A. In our observa¬ 
tions we normally split the total bandwidth into 20 parts and then 
process them separately, each with their own polyco files. This 
helps to mitigate the effect, giving a maximum TOA error for the 
lowest-frequency part on the order of ~ 10 ps for most MSPs, 
with the most problematic pulsars being J0737-3039A (0.64 ms 
residual smearing), J2215-1-5135 (0.16 ms), J1816-I-4510 (55 ps), 
J0218-H4232 (52 ps), J0621-H1002 (46 ps), and J0751-H1807 
(39 ps). The largest error in terms of spin phase occurs for PSR 
J2215+5135, where the effect is equivalent to 0.06 rotational cy¬ 
cles. Eor most of our observed MSPs, this error is about 0.01 
cycles or less. 

3.1.2. Double pulsar J0737-3039A 

Despite its low declination and maximum elevation of 6.4° (as 
seen from the LOEAR Core), we detected the 22-ms A-pulsar 
in the double pulsar system J0737-3039 in a single 1-h obser¬ 
vation (see Eig. 1). To exclude any broadening due to the ef¬ 
fect described in Sect. 3.1.1 we have also reprocessed the data, 
folding each channel separately. It did not make any noticeable 
difference to the profile, because at the epoch of observation the 
pulsar had a relatively low orbital velocity (orbital phase ~0-0.4, 
longitude of periastron ~ 300.5°). Erom high-frequency observa¬ 
tions it is known that the B-pulsar is not visible anymore due to 
large orbital (17 °/yr; Kramer et al. 2006a) and geodetic (5 °/yr; 
Breton et al. 2008) precession, but is expected to reappear again 
in ~ 2035 (Perera et al. 2010). However, we might be able to 
catch its emission at our low frequencies with LOEAR, if a wider 
low-frequency beam (Cordes 1978b) is still grazing our line-of- 
sight. Thus, we folded our data with the B-pulsar ephemeris as 
well, but we did not detect it after also searching in period and 
DM. Radio emission from the B-pulsar was only ever visible 
at two narrow orbital phase windows (~ 0.58 and ~ 0.78) when 
emission from the A-pulsar affected the magnetosphere of the B- 

http://tempo.sourceforge.net/ 
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Table 2. List of the 48 detected MSPs. 


PSR 

P 

(ms) 

DM 

(pc cm“^) 

Binary? 

WSRT 

detected? 

BSA 

detected? 

LOFAR 

ObsID 

Epoch 

(MJD) 

Fobs 

(min) 

Prof. 

Sign. 

S/N* 

LBA 

detected? 

Ref 

J0030-H0451 

4.865 

4.333 

Isolated 


y 

L83021 

56304.694 

20 

251 

31 

y 

1 

J0034-0534 

1.877 

13.765 

He WD 

y 

y 

L81272 

56286.738 

20 

543 

63 

y 

1 

J0214-H5222 

24.575 

22.037 

He WD/sdB? 



L196378 

56646.791 

20^' 

116' 

411' 


4 

J0218-H4232 

2.323 

61.252 

He WD 

n 

y 

L155442 

56473.304 

20 

126 

19 


1 

J0337-H1715 

2.733 

21.316 

He WD-rHe WD 



L167133 

56512.229 

60 

50 

11 


2 

J0407-H1607 

25.702 

35.65 

He WD 



L227494 

56790.515 

20 

194 

30 


20 

J0621-H1002 

28.854 

36.601 

CO WD 

n 

y 

L81270 

56289.023 

20 

47 

9 


3 

J0636+5129 

2.869 

11.107 

BW? 



L196371 

56648.046 

lO'^ 

301' 

lOi' 


4 

J0645+5158 

8.853 

18.247 

Isolated 



L85909 

56322.009 

20 

59 

19 


4 

J0737-3039A 

22.699 

48.920 

PSR 



L85911 

56321.940 

60 

25 

6 


5 

J0751-H1807 

3.479 

30.249 

He WD 


y 

L81051 

56280.049 

20 

24 

7 


1 

J1012-H5307 

5.256 

9.023 

He WD 

y 

y 

L81268 

56289.149 

20 

183 

32 

n 

1 

J1022-H1001 

16.453 

10.252 

CO WD 

y 

y 

L81254 

56296.126 

20 

274 

49 

n 

6 

J1023-H0038 

1.688 

14.325 

MS/Redback 

y 


L85233 

56315.163 

20 

105 

30 


7 

J1024-0719 

5.162 

6.485 

Isolated 

n 

y 

L81049 

56280.168 

20 

40 

10 

n 

1 

J1038+0032 

28.852 

26.59 

Isolated 



L227490 

56782.804 

20 

39 

8 


21 

J1231-1411 

3.684 

8.090 

He WD 



L227492 

56789.862 

20 

27 

7 


17 

B1257H-12 

6.219 

10.166 

Planets 

y 

y 

L81253 

56296.230 

20 

313 

44 

n 

8 

J1453-H1902 

5.792 

14.049 

Isolated 



L227293 

56779.988 

20 

18 

6 


22 

J1544-H4937 

2.159 

23.226 

BW 



L227294 

56780.003 

20 

89 

23 


23 

J1640-H2224 

3.163 

18.426 

He WD 


y 

L81266 

56289.368 

20 

66 

15 


1 

J1709-H2313 

4.631 

25.347 

He WD 



L249810 

56964.533 

20 

14 

5 


24 

J1713-H0747 

4.570 

15.992 

He WD 

n 

y 

L149156 

56465.941 

60 

30 

8 


1 

J1730-2304 

8.123 

9.617 

Isolated 


y 

L164998 

56490.889 

30 

32 

9 


1 


Notes. P and DM are from the best ephemerides we used. Col. 4 shows the binary companion or system type, where ‘BW’ denotes ‘black widow’ eclipsing systems with ultra-light (UL) companion 
(M < O.OSMq), ‘Redback’ - compact eclipsing systems with more massive companion (M ~ 0.1-0.2Mq), ‘MS’ - main-sequence star companion, ‘WD’ - white dwarf companion, either helium 
(He), or CO/ONeMg (CO), ‘sdB’ - subdwarf B star. WSRT observations are from Stappers et al. (2008) and from Archibald et al. (2009) in the case of PSR J1023-I-0038. BSA is the Large Phased 
Array in Pushchino and observations are those reported by Kuzmin & Losovsky (2001). The corresponding flux density measurements are given in Table 4 for those pulsars detected with LOFAR. 
The LOFAR ObsID, Epoch, Tabs, S/N and profile significance (Prof. Sign.) are for the best individual observation. Profile significance is defined as in Lorimer & Kramer (2005, p. 167, Eq. 7.1). 
The last column gives a reference to the ephemerides used. P' Corresponding value is for the sum of several 5-min observations. **) Given the number of profile bins as in Eig. 1. 

References. (1) Desvignes & et al. (2015); (2) Ransom et al. (2014); (3) Splaver et al. (2002); (4) Stovall et al. (2014a); (5) M. Kramer (private communication); (6) Verbiest et al. (2009); 
(7) Archibald et al. (2009); (8) Konacki & Wolszczan (2003); (9) J. Hessels (private communication); (10) Gonzalez et al. (2011); (11) Toscano et al. (1999); (12) Lynch et al. (2013); (13) Champion 
et al. (2005); (14) Arzoumanian et al. (1994); (15) Guillemot et al. (2012); (16) Doroshenko et al. (2001); (17) Ransom et al. (2011); (18) Camilo et al. (1996b); (19) Cognard et al. (2011); 
(20) Lorimer et al. (2005); (21) Burgay et al. (2006); (22) Lorimer et al. (2007); (23) Bhattacharyya et al. (2013); (24) Lewandowski et al. (2004); (25) Kaplan et al. (2013). 
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Table 2. continued. 


PSR 

P 

(ms) 

DM 

(pc cm“^) 

Binary? 

WSRT 

detected? 

BSA 

detected? 

LOFAR 

ObsID 

Epoch 

(MJD) 

Lobs 

(min) 

Prof. 

Sign. 

S/N* 

LBA 

detected? 

Ref 

J1738-H0333 

5.850 

33.778 

He WD 



L124889 

56399.183 

20 

21 

6 


1 

J1744-1134 

4.075 

3.139 

Isolated 

y 

y 

L81264 

56293.440 

20 

37 

15 


1 

J1810-H1744 

1.663 

39.659 

BW 



L81263 

56293.456 

20 

861 

65 

n 

9 

J1816-H4510 

3.193 

38.887 

He? WD/sdB? 



L212736 

56737.245 

20'' 

8881' 

159+ 


4, 25 

J1853-H1303 

4.092 

30.570 

He WD 



L84523 

56311.438 

20 

24 

7 


10 

B18554-09 

5.362 

13.300 

He WD 


y 

L131365 

56417.151 

20 

23 

5 


1 

J1905-H0400 

3.784 

25.692 

Isolated 



L249826 

56964.657 

20 

18 

7 


10 

J1911-1114 

3.626 

30.975 

He WD 

y 

y 

L81277 

56286.527 

20 

25 

8 


11 

J1918-0642 

7.646 

26.554 

He WD 



L81276 

56286.543 

20 

28 

7 


1 

J1923-H2515 

3.788 

18.858 

Isolated 



L85594 

56318.472 

20 

58 

12 


12 

B1937H-21 

1.558 

71.040 

Isolated 



LI 38647 

56434.134 

30 

275 

25 


1 

J1944-H0907 

5.185 

24.34 

Isolated 



L84521 

56311.472 

20 

103 

12 


13 

B1953H-29 

6.133 

104.501 

He WD 



L84522 

56311.456 

20 

44 

6 


10 

B19574-20 

1.607 

29.117 

BW 

y 


L81275 

56286.559 

20 

230 

24 


14 

J2019-H2425 

3.935 

17.203 

He WD 


y 

L146225 

56457.103 

20 

18 

6 


1 

J2043-H1711 

2.380 

20.710 

He WD 



L84518 

56311.522 

20 

96 

14 


15 

J2051-0827 

4.509 

20.745 

BW 

n 

y 

L85592 

56318.504 

20 

65 

12 


16 

J2145-0750 

16.052 

8.998 

CO WD 

y 

y 

L81259 

56293.607 

20 

324 

37 

y 

1 

J2214-H3000 

3.119 

22.557 

BW 



L146228 

56457.159 

20 

27 

7 


17 

J2215-H5135 

2.610 

69.2 

MS/Redback 



L85588 

56318.567 

20 

114 

13 


9 

J2235-H1506 

59.767 

18.09 

Isolated 


y 

L168068 

56521.018 

30 

25 

9 


18 

J2302-H4442 

5.192 

13.762 

He(?) WD 



L84516 

56311.603 

20 

40 

10 


19 

J2317-H1439 

3.445 

21.907 

He WD 


y 

L83022 

56304.635 

20 

176 

40 

n 

1 

J2322-H2057 

4.808 

13.372 

Isolated 


y 

L146234 

56460.218 

20 

14 

9 


1 


Notes. P and DM are from the best ephemerides we used. Col. 4 shows the binary companion or system type, where ‘BW’ denotes ‘black widow’ eclipsing systems with ultra-light (UL) companion 
(M < O.OSMq), ‘Redback’ - compact eclipsing systems with more massive companion (M ~ O.I-O. 2 M 0 ), ‘MS’ - main-sequence star companion, ‘WD’ - white dwarf companion, either helium 
(He), or CO/ONeMg (CO), ‘sdB’ - subdwarf B star. WSRT observations are from Stappers et al. (2008) and from Archibald et al. (2009) in the case of PSR J1023-I-0038. BSA is the Large Phased 
Array in Pushchino and observations are those reported by Kuzmin & Losovsky (2001). The corresponding flux density measurements are given in Table 4 for those pulsars detected with LOFAR. 
The LOFAR ObsID, Epoch, Tobs. S/N and profile significance (Prof. Sign.) are for the best individual observation. Profile significance is defined as in Lorimer & Kramer (2005, p. 167, Eq. 7.1). 
The lasf column gives a reference to the ephemerides used. Corresponding value is for the sum of several 5-min observations. Given the number of profile bins as in Eig. 1. 

References. (1) Desvignes & et al. (2015); (2) Ransom et al. (2014); (3) Splaver et al. (2002); (4) Stovall et al. (2014a); (5) M. Kramer (private communication); (6) Verbiest et al. (2009); 
(7) Archibald et al. (2009); (8) Konacki & Wolszczan (2003); (9) J. Hessels (private communication); (10) Gonzalez et al. (2011); (11) Toscano et al. (1999); (12) Lynch et al. (2013); (13) Champion 
et al. (2005); (14) Arzoumanian et al. (1994); (15) Guillemot et al. (2012); (16) Doroshenko et al. (2001); (17) Ransom et al. (2011); (18) Camilo et al. (1996b); (19) Cognard et al. (2011); 
(20) Lorimer et al. (2005); (21) Burgay et al. (2006); (22) Lorimer et al. (2007); (23) Bhattacharyya et al. (2013); (24) Lewandowski et al. (2004); (25) Kaplan et al. (2013). 
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Table 3. List of the 27 non-detected MSPs. 


PSR 

P 

(ms) 

DM 

(pc cm“^) 

Binary? 

WSRT 

detected? 

BSA 

detected? 

LOFAR 

ObsID 

Epoch 

(MJD) 

Fobs 

(min) 

^ Urn 

(mJy) 

KLOl flux density 
(mJy @ 102 MHz) 

Ref 

J0023+0923 

3.050 

14.326 

BW 



L84514 

56311.679 

20 

4 


1 







LI 55443 

56473.224 

20 




J0340+4130 

3.299 

49.575 

Isolated 



L85586 

56318.785 

20 

3 


9 







LI 55441 

56473.360 

20 




J0613-0200 

3.062 

38.781 

He WD 

n 

y 

L81271 

56289.007 

20 

5 

240 + 100 

1 







LI 55440 

56473.467 

20 




J1327-0755 

2.678 

27.912 

He WD 



L243060 

56907.580 

20 

5 


29 

J1400-1438* 

3.084 

4.928 

He WD 



L173398 

56533.626 

20 

10 


42 

J1455-3330 

7.987 

13.562 

He WD 



L85908 

56322.243 

20 

60 


1 

J1551-0658 

7.09 

21.6 

Binary? 



L227493 

56789.992 

20 

7 


30 

J1600-3053 

3.598 

52.317 

He WD 



L146215 

56456.919 

20 

70 


1 

J1614-2230 

3.151 

34.487 

CO WD 



L85907 

56322.269 

20 

25 


25 

B1620-26 

11.076 

62.865 

He WD+UL 


y 

L85906 

56322.285 

20 

40 

450 + 200^ 

1 

J1643-1224 

4.622 

62.377 

He WD 



L81278 

56286.424 

20 

15 


1 







L85231 

56315.324 

20 




J1646-2142 

5.85 

29.8 

Binary 



L227491 

56787.043 

20 

30 


31 

J1719-1438 

5.79 

36.862 

UL/Carbon WD? 



L249812 

56964.547 

20 

15 


32, 33 

J1721-2457 

3.497 

47.758 

Isolated 



L249814 

56964.562 

20 

45 


34, 35 

J1741+1351 

3.747 

24.0 

He WD 



L81281 

56281.492 

20 

4 


26 







L83025 

56304.456 

20 




J1745+1017 

2.652 

23.970 

BW? 



L249816 

56964.576 

20 

4 


36 

J1751-2857 

3.915 

42.808 

He WD 



L249818 

56964.591 

20 

220 


37 

J1804-2717 

9.343 

24.674 

He WD 



L249820 

56964.606 

20 

no 


38, 39 

J1828+0625 

3.63 

22.4 

Binary 



L249822 

56964.620 

20 

6 


31 

J1832-0836 

2.719 

28.18 

Isolated 



L249824 

56964.635 

20 

40 


40 

J1903+0327 

2.150 

297.548 

MS 



L85596 

56318.440 

20 

20 


1 

J1910+1256 

4.984 

38.097 

He WD 



L85595 

56318.456 

20 

10 


1 

J1911+1347 

4.626 

30.99 

Isolated 



L243072 

56907.790 

20 

5 


41 

J1949+3106 

13.138 

164.126 

CO WD 



L85593 

56318.487 

20 

4 


27 

J2010-1323 

5.223 

22.181 

Isolated 



L84520 

56311.490 

20 

15 


1 

J2017+0603 

2.896 

23.918 

He WD 



L85591 

56318.520 

20 

5 


19 

J2229+2643 

2.978 

23.019 

He WD 


y 

L84517 

56311.587 

20 

3 

60 + 30 

28 


Notes. Column descriptions are the same as in Table 2. S um is the upper limit on the mean flux density and taken as three times the nominal error of the mean flux density measurements. Flux 
densities hy Kuzmin & Losovsky (2001, KLOl) are at 102 or 111 MHz. Some of the pulsars were attempted twice. Observations were done both at 102 and 111 MHz, and it is not clear for which 
frequency the flux densities were reported. 1*1 The catalogued position of PSR J1400-1438 was about 6' off from the position derived from over three years of timing with the GET (J. Swiggum, 
private communication). We detected this MSP in a recent LOFAR observation using the updated coordinates. The average profile, effective pulse width and flux density measurement for this MSP 
will be presented elsewhere. 

References. (1) Desvignes & et al. (2015); (9) J. Hessels (private communication); (19) Cognard et al. (2011); (25) Demorest et al. (2010); (26) Jacoby et al. (2007); (27) Deneva et al. (2012); 
(28) Wolszczan et al. (2000); (29) Boyles et al. (2013); (30) Hessels et al. (2011); (31) Ray et al. (2012); (32) Ng et al. (2014); (33) Bailes et al. (2011); (34) Edwards & Bailes (2001); (35) Janssen 
et al. (2010); (36) Barr et al. (2013b); (37) Stairs et al. (2005); (38) Hobbs et al. (2004b); (39) Lorimer et al. (1996); (40) Burgay et al. (2013); (41) Lorimer et al. (2006); (42) Rosen et al. (2013). 
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J0645+5158 Nbm=512 





0.2 0,4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1 


J1023+0038 I Nbln=12B . J1024-0719 , Nbm=12B 





0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 O.B 1 0 0.2 0.4 0.6 O.B 1 0 0.2 0.4 






0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 O.B 1 C 0.2 0.4 0.6 0.8 1 

J1944+0907 L Nbln=256 



0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 O.B 1 0 0.2 0.4 0.6 O.B 1 0 0.2 0.4 0.6 O.B 1 







0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1 

Fig. 1. LOFAR total-intensity pulse profiles for all detected MSPs from the best (highest S/N) individual observations. The horizontal axes show 
the pulse phase, and the vertical axes are the flux density in arbitrary units. 
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pulsar (Lyne et al. 2004), and we missed these windows in our 
observation. Therefore, more observations are needed to rule out 
emission from the B-pulsar at low frequencies. 

3.1.3. PSRJ1023+0038 

We detected the 1.69-ms pulsar J1023+0038, which is in a 4.7-h 
orbit with a bloated stellar companion and has been called the 
‘missing link’ between low-mass X-ray binaries (LMXBs) and 
radio MSPs (Archibald et al. 2009). We observed it five times, 
with 20-min integration times per observation, before February 
2013 and detected it twice in the orbital phase ranges 0.64-0.71 
and 0.84-0.91. The pulsar was very likely eclipsed in the other 
three observations (at the observed frequencies, it is expected 
to be eclipsed for at least orbital phases 0.0-0.5, see Archibald 
et al. 2013). The profile from the best observation on January 23, 
2013 is shown in Fig. 1. The profile has a two-component shape, 
as seen in WSRT observations at 155-157 MHz by Archibald 
et al. (2009). Two observations on August 31, 2013 and October 
5, 2013 at non-eclipse orbital phases did not show any emis¬ 
sion from PSR J1023+0038. This was confirmed earlier with 
high-frequency observations using other radio telescopes. This 
is evidence that it switched again into its accreting LMXB phase 
(Patruno et al. 2014; Stappers et al. 2014). 

3.1.4. PSR B1937+21 

We successfully detected the first millisecond pulsar, B1937+21, 
in the LOFAR HBA band. This is the first ever published detec¬ 
tion of the pulse profile of this pulsar at these low frequencies 
(see Fig. 2). The main pulse (MP) and interpulse (IP) compo¬ 
nents are clearly distinguished at the highest frequencies. Fig¬ 
ure 2 shows the profile evolution in the HBA band from 105 to 
197 MHz. The rapid change in the profile shape due to scattering 
is evident. Below 150 MHz the IP is no longer visible and below 
120 MHz the whole profile is completely washed out by scatter¬ 
ing. As mentioned before, scattering time, Tscat, scales as 
with frequency, /. Following the empirical relationship between 
Tscat, DM and / from Bhat et al. (2004), the Tjcat should increase 
from ~ 2 to ~ 19 ms when going from 197 to 105 MHz, i.e. from 
~ 1 to almost ~ 12 spin periods {P - 1.56 ms). Joshi & Kramer 
(2009) reported on the IP giant-pulse component in the average 
profile from their observations at a frequency of 238 MHz with 
the Giant Metrewave Radio Telescope (GMRT). We do not see 
this giant-pulse-related hump in the average profile of either the 
MP or the IP components, most likely because of much larger 
profile broadening. Kuzmin & Losovsky (2002) reported the de¬ 
tection of four individual giant pulses at 112 MHz, though the 
pulse trains they present are very noisy. We have also performed 
a search for strong giant pulses in the full band, but did not find 
anything significant. Only searching the high-frequency half of 
the band or using the cyclic spectroscopy technique (Demorest 
2011; Walker et al. 2013; Archibald et al. 2014) may be needed. 
We defer more detailed single-pulse study of PSR B1937+21 for 
a subsequent paper. 

3.2. Flux densities 
3.2.1. Contributing factors 

For a number of reasons, flux calibration of beamformed LO¬ 
FAR data is non-trivial. The main uncertainty comes from the 
fact that the characteristic beam shape and sensitivity strongly 
depend on the elevation and azimuth of the source. LOFAR also 



Pulse phase 


Fig. 2. Profile evolution of pulsar B1937+21 within the LOFAR HBA 
band in the range 105-197 MHz. Each profile has 128 bins and is folded 
in ~ 9.2-MHz wide sub-bands centred at the labelled frequency in MHz. 

has a very large fractional bandwidth, thus the frequency depen¬ 
dence of the beam size and differences in the sky temperature 
across the bandpass should be considered. Overall, the follow¬ 
ing aspects have to be taken into account: 

Effective area. The total effective collecting area of the LO¬ 
FAR array, Aefr, or gain, G, is frequency dependent. For a 48- 
tile HBA station the upper limit on the effective area, a™’'(/I), 
is given by van Haarlem et al. (2013)'^: 

a”ff“(4) = min {2564^, 1200} m^, (1) 

where A is the observing wavelength in m. 

Beam model. An accurate beam model of the LOFAR HBAs at 
110-190 MHz provides the dependence of the antenna gain, 

''' http: //WWW. astron.nl/radio-observatory/astronomers/ 
lofar-imaging-capabilities-sensitivity/sensitivity- 
lofar-array/sensiti 
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G, on the zenith angle and azimuth for a single tile or station, 
and preferably takes into account the shadowing of individ¬ 
ual tiles within a station at low elevations. 

System temperature. The system temperature, Tsys = Tsky + 
Ta, is also frequency dependent. The main contribution to 
Tsys comes from the sky temperature, Tsky, which strongly 
changes with frequency as (Lawson et al. 1987). We 
should note, that, although the spectral index is remarkably 
uniform for most of the Galactic emission, it does vary across 
the sky in the range from -2.40 to -2.55 (Roger et al. 1999). 
The antenna temperature, Ta, is not constant either. At the 
centre of the band its average value is about 400 K (Arts et al. 
2013), although it increases slightly towards higher frequen¬ 
cies (van Haarlem et al. 2013). Below 120MHz and above 
180 MHz Ta increases rapidly. We used a sixth-order poly¬ 
nomial function of frequency based on Tsys measurements 
from Wijnholds & van Cappellen (2011) to fit for Ta and 
calculate its value for each frequency channel. 

Broken tiles. In practice, not all the HBA tiles are operational 
at all times, and on average about 5-10% of them are unused 
while observing, making the effective area smaller. For some 
of the earlier observations this fraction was as high as 30%. 

Coherent station summation. Coherent summation of the sta¬ 
tions should increase the sensitivity linearly with the number 
of stations added. In practise, however, we found the depen¬ 
dence of the S/N to scale as Nj, where W is the number of 
48-tile stations, and y = 0.85 is the coherence factor. The 
possible reasons for this non-linear dependence are corre¬ 
lated noise between stations and uncalibrated non-zero phase 
offsets in the phase-frequency dependence for each of the 
polarisations/stations, which can be caused by the station’s 
electronics themselves or the influence of the ionosphere. 
This measurement was based on the BG/P data, and we used 
BG/P for most of our observations. Similar tests are currently 
being done for the Cobalt correlator, where these non-zero 
phase offsets will be taken into account. 

Radio frequency interference. RFI clipping reduces the total 
on-source time and bandwidth. On average, 25-30% of our 
data were clipped, the smallest and largest fractions being 
10% and 46%, respectively. This large zapping fraction is an 
unfortunate consequence of keeping a frequency resolution 
of only 195 kHz (this is needed in order to achieve the desired 
time resolution). The typical spectral width of the RFI seen 
by LOFAR is ~ 1 kHz, and thus higher-frequency-resolution 
data can be excised of RFI without sacrificing such a large 
fraction of the data. 

3.2.2. Calibration 

The Stokes vector. So = {Iq, Qq, Uq, Vq), of the incident radia¬ 
tion is related to the measured Stokes vector, S = (/, Q, U, V), 
via a 4 X 4 Mueller matrix, M, as S = M ■ Sq. The measured 
Stokes I component will depend on the initial Stokes So as 
I - niiilo -H nijQQo + mjijUo -I- mjvVo, where mjj are correspond¬ 
ing elements of M. In our observations mjj is the dominant fac¬ 
tor and we can assume that all observed power comes in Stokes 
I. Pulsars are predominantly very weakly circularly polarised, 
thus nijv is small compared to mu (< 1%). For linear polarisa¬ 
tion, niiQ and niju are of the order 5-10% and strongly depend 
on the hour angle (Noutsos et al. 2015). Noutsos et al. (2015) 
also report significant leakage to circular polarisation at large 
hour angles, however all our observations were carried out near 
transit where the leakage is minimal. Moreover, Faraday rotation 
from the ISM removes most of the contribution from linear po¬ 


larisation in our large observing band. Therefore, we can simply 
consider I » mijlo. 

In our flux density measurements we used the Hamaker beam 
model (Hamaker 2006, and references therein), which is also 
used in the LOFAR software for initial phase calibration of in¬ 
terferometric imaging data (van Haarlem et al. 2013). Based on 
this model we calculated the Jones matrices of antenna response 
for a given HBA station, frequency and sky direction using the 
mscorpol package'^. Then, Jones matrices (e.g. Jones 1941; Fy- 
mat 1971) can be transformed to a corresponding Mueller ma¬ 
trix, and in particular for the mu'. 

niuif ,Z,A) - — [JxxJxx -^xyJxij + JyxJyx -^yy-^yy) ’ (2) 

where Jij are components of the Jones matrices, / is the fre¬ 
quency, z is the zenith angle, and A is the azimuth (e.g. Hamaker 
et al. 1996). Finally, to link with the absolute flux scale we nor¬ 
malised our mu(f,z,A) by the corresponding values for a refer¬ 
ence observation of Cassiopeia A by Wijnholds & van Cappellen 
(2011). In calculations of Jones matrices we always used the sta¬ 
tion CSOOl as our template. The difference in flux densities when 
other stations are used to derive Jones matrices is much smaller 
than the nominal flux uncertainty. Using niu, the fleft for a given 
frequency and direction can be simply calculated as: 

«eff(/, z, A) = X '«//(/> A). (3) 

To determine the pulse-phase-averaged flux density of the 
pulsar (hereafter referred to as ‘mean flux density’), we adapt 
Equation (7.10) of Lorimer & Kramer (2005) to obtain the flux 
density, S ;, of the /-th profile bin: 

Si = (S/N),- ^ (4) 

‘-f V”p(7’ote/A'b)B 

where /? is a digitisation correction factor, (S/N), is the signal- 
to-noise ratio of the /-th profile bin, Up = 2 is the number of 
polarisations summed. Tabs is the observation length, B is the 
observation bandwidth, Nb is the number of profile bins, G = 
Agffllk is the telescope gain, Aeff is the total effective area, k is 
the Boltzmann constant, and Tsys is the system temperature. For 
systems with a large number of bits per sample (8 in our case) 

Expanding this expression to apply to LOEAR, for the flux 
density S i{f, z. A) in Jy of the /-th profile bin of a given frequency 
channel and sub-integration, we derive: 

^ 2;6k(S/N)f^-^[7’A(/) + r.ky)/, /, b)] 

^ Nla,ff{f,z,A)muif,z,A)[l-^] 


^n,[l - am^)Af' 

where I and b are Galactic longitude and latitude, ^ is the overall 
fraction of bad tiles, ^(/) is the RFI fraction. Tint is the length 
of a sub-integration. A/ is the width of the frequency channel, 
and (S/N)J’^’^ = (JfpA _ (X))/crx is the signal-to-noise ratio 
of the /-th profile bin of the individual frequency channel and 
sub-integration, with Xz ’ the value of the /-th profile bin, {X} 
the mean and ctx the rms value (and both {X) and ctx are cal¬ 
culated in the off-pulse window). To calculate Tsky we use the 
reference sky temperature at 408 MHz from the skymap deter¬ 
mined by Haslam et al. (1982) and scale it to the frequency of 

https://github.com/2baOrNot2ba/mscorpol 
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the channel using (Lawson et al. 1987). At present there 
are few sophisticated sky models available (de Oliveira-Costa 
et al. 2008; Polisensky 2007), which are based on published ref¬ 
erence sky maps, and that we can use in the future especially for 
the calibration in the LBA frequency range 10-90 MHz (see also 
Nelles & et al. 2015). 

Finally, to measure the mean flux density 5 mean we find the 
average value of S A) as: 

, Nb M 

^™ean=^EE Si(f{j),z(j),A{j)), (6) 

/=! M 

where M = A^chA)nt - A^zap with Ach the total number of fre¬ 
quency channels. Aim the total number of sub-integrations, and 
Azap the number of completely zapped frequency channels and 
sub-integrations. 

We estimated the nominal uncertainty, cr,(/, z,A), on 

f z A 

S i{f,z,A) by taking (S/N),’ ’ = 1 for every individual channel 
and sub-integration, i.e. crj{f,z,A) - Si(f,z,A)\^^^j^^t.i.A_j^. Then, 
we And the nominal uncertainty, iTmean, on the mean flux density, 

5mean5 


M^M 


,crHf(j),z(j),A(j)) 


(7) 


However, there are also other factors that can affect our flux 
measurements, but which are very difficult (or not possible) to 
take account of, and were not addressed here. For completeness 
we note that these are: 


Scattering. Our flux density measurements can be seriously af¬ 
fected if Tscat is close to or exceeds a pulse period. In these 
cases, the measured S/N is reduced and we underestimate 
Tsys due to the contribution from the scattered pulse signal. 

Refractive scintillation. Diffractive scintillation should not af¬ 
fect our measurements because the decorrelation bandwidths 
are lower than the width of a channel (195 kHz) for all of 
our sources, so many scintles are averaged out. On the other 
hand, refractive interstellar scintillation (RISS) can change 
the pulsar flux density by a factor of ~ 1.5 at low observ¬ 
ing frequencies. Gupta et al. (1993) found the characteristic 
RISS time scale to vary from a few days to several weeks 
at a frequency of 74 MHz for a sample of nearby pulsars. 
The measured modulation indices (ratio of the standard de¬ 
viation of the observed flux densities to their mean) ranged 
from 0.15-0.45. The RISS time scale and modulation index 
scale with the observing frequency as oc and oc 
respectively (Rickett et al. 1984; Coles et al. 1987). Obser¬ 
vations spanning dozens of months are required to measure 
the average flux densities more accurately. Here we provide 
only the flux density measurements for the best individual 
observation for each pulsar. 

Beam jitter. The jitter of the tied-array beam is caused by 
the ionosphere and can be as large as ~ 2-3', equal to the 
half-width at half-maximum of the HBA full-core tied-array 
beam^®, but is usually much smaller than this. However, 
when the ionosphere is excited, e.g. around sunset/sunrise, 
it could move the source to the edge of the beam leading to 
inferred flux densities lower by a factor of ~ 2. In future anal¬ 
yses forming a ring of tied-array beams around the pointing 

For a demonstration, see: http://www. astron.nl/dailyimage/ 
index.html?main.php?date=2014012 3 


position could eliminate this problem by choosing the beam 
with highest S/N in the post-processing. This, however, is 
expensive with respect to both data volume and processing 
time. 

Galactic plane. Variation of the system temperature with time 
due to rise/set of the Galactic plane can increase Tsys by 30- 
40% if it is in the field-of-view (FoV) of the primary sta¬ 
tion beam or grating lobes. This is also possible for other 
strong background sources in the FoV. Also, Tsys can vary 
with pointing direction due to noise coupling effects. 

To get realistic values for the systematic uncertainty on our 
flux density measurements, we also estimated the flux density 
using images from the Multifrequency Snapshot Sky Survey 
(MSSS, Heald et al. 2015) for a sample of bright MSPs and slow 
pulsars that we could unambiguously identify. On average our 
mean flux density measurements agree with those from MSSS 
data within ~ 40%. The same calibration technique was used 
to get several monthly flux density measurements for ten bright 
non-millisecond pulsars in the LOFAR HBA range (Bilous et 
al. 2015, in prep.). Comparing measured fluxes to spectra fit 
through the literature values, the authors estimate the uncertainty 
of a single LOFAR flux measurement as 50%. Thus, MSP fluxes 
reported here should also be used with this uncertainty. 

3.2.3. Flux measurements 

Table 4 lists the measured 110-188 MHz flux densities of our 
MSP sample. In the Table we give the catalogued value of DM 
for reference, followed by our measurements of effective pulse 
width, Weff (Col. 3), duty cycle, 6 = WesIP (Col. 4), and mean 
flux density measurements, 5mean (Col. 5). Given 5 mean errors in 
parentheses are nominal uncertainties only from the Eq. 7. For 
each MSP the (S/N)J’^’^ was calculated based on (X) and ctx 
in the pre-selected off-pulse window. For MSPs with broad pro¬ 
files (due to scattering) we tried to place the off-pulse window 
just before the leading edge of the profile, but most likely this 
still resulted in underestimated S/N and flux density. The Wefr 
was calculated as the integrated flux density in the pulse pro¬ 
file over its peak flux density. In Col. 6 we list the flux density 
measurements at 102/111 MHz reported by Kuzmin & Losovsky 
(2001). In Col. 7 we provide the predicted mean flux density us¬ 
ing the reference flux density value, 5 ref (Col. 8) at the reference 
observing frequency, Vref (Col. 9), and spectral index (Col. 10) 
from published data. We give references in Col. 11 to all the 
published data we used in our calculations of the predicted flux 
density. The values of the predicted mean flux densities listed 
are the average between the values at the edges of our frequency 
range of 110-188 MHz. If the spectral index was not known, we 
used the mean value of -1.4 + 1 for the pulsar population (Bates 
et al. 2013). We used this recently published value instead of the 
generally accepted value of -1.8 + 0.2 by Maron et al. (2000) as 
it takes into account the sample bias in the spectral index mea¬ 
surements available. Although the work of B ates et al. (2013) fo¬ 
cused on the normal pulsar population, MSP spectra seem to be 
similar to those of normal pulsars (Kramer et al. 1998), and thus 
we follow that assumption here. For the non-detected MSPs, we 
provide upper limits 5iim (Table 3). We calculated 5iim as three 
times the nominal error on the mean flux density measurement. 

Figure 3 shows the comparison between the measured and 
predicted flux densities for our 48 detected MSPs. Predicted 
flux densities were calculated from the available high-frequency 
data and known spectral indices as given in Table 4. For at 
least a third of the MSPs the main uncertainty on the pre- 
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Table 4. Mean flux densities for the detected MSPs in the frequency range 110-188 MHz. 


PSR 

DM 

(pc cm“^) 

IVefF 

(ms) 

b 

(%) 

^ mean 

(mJy) 

KLOl flux density 
(mJy @ 102 MHz) 

Predicted S mean 

(mJy) 

(mJy) 

Vref 

(MHz) 

Spectral 

index 

Reference 

J0030+0451 

4.333 

0.61(1)* 

12.6(2)* 

99(1) 

380 + 200'' 

80 + 40 

7.9+ 0.2 

430 

-2 1 

1, 38 

J0034-0534 

13.765 

0.54(1)* 

28.8(4)* 

491(2) 

250+ 120 

400 + 100 

17 + 5 

436 

-2.64 + 0.05 

2, 38 

J0214+5222 

22.037 

1.5(2)* 

6.1(8)* 

21.0(7) 


4.1 

0.90 

350 

-1.6 

32 

J0218+4232 

61.252 

0.77(2) 

33.1(8) 

37.6(8) 

270+ 150* 

530 + 70 

30^0 

410 

-2.41 + 0.04 

3, 38 

J0337+1715 

21.316 

0.85(4)* 

16(2) 

5.0(7) 


50 + 100 

~2 

1400 

-1.4 + 1 

30, 22 

J0407+1607 

35.65 

4.2(1)* 

16.4(4)* 

56.5(8) 


50 + 50 

10.2 

430 

-1.4 + 1 

34, 22 

J0621+1002 

36.601 

6.1(2) 

21.3(8) 

20(1) 

50 + 25'' 

170 + 380 

1.9+ 0.3 

800 

-2.5 + 1.4 

4 

J0636+5129 

11.107 

0.14(2) 

4.9(8) 

4.1(8) 


5.5 

1.8 

350 

-1.2 

32 

J0645+5158 

18.247 

0.16(2)* 

1.8(2)* 

4.7(7) 


20 

2.4 

350 

-2.2 

32 

J0737-3039A 

48.920 

2.5(2) 

11.2(8) 

64(10) 


40 + 80 

1.6+ 0.3 

1400 

-1.4 + 1 

5,22 

J0751+1807 

30.249 

0.65(5) 

19(2) 

5.5(8) 

70 + 30 

30+10 

10 

430 

-0.9 +0.3 

6, 4 

J1012+5307 

9.023 

0.69(2)* 

13.1(4)* 

34.5(5) 

30+15 

220 + 90 

30 

410 

-1.8 +0.4 

7, 4 

J1022+1001 

10.252 

1.00(3)* 

6.1(2)* 

40.7(9) 

90 + 40 

110 + 40 

20 + 9 

370 

-1.7 +0.1 

8, 4 

J1023+0038 

14.325 

0.15(1)* 

9.1(8)* 

48(2) 


440-3400 

0.4-3.1*^ 

1600 

—2.8 

31, 9 

J1024-0719 

6.485 

0.54(4)* 

10.4(8)* 

21(2) 

200+ 100 

30+10 

4.6+ 1.2 

410/430 

-1.5 +0.2 

10, 4 

J1038+0032 

26.59 

5.4(2) 

18.7(8) 

16(1) 


5 + 10 

0.20 + 0.04 

1400 

-1.4 + 1 

5,22 

J1231-1411 

8.090 

0.41(3) 

11.1(8) 

22(3) 


5 + 7 

0.4 

820 

-1.4 + 1 

27, 22 

B1257+12 

10.166 

0.63(1)* 

10.2(2)* 

80.9(9) 

150 + 50 

180 + 90 

20 

430 

-1.9 +0.5 

11, 4 

J1453+1902 

14.049 

0.96(9) 

17(2) 

6(1) 


10+10 

2.2+ 0.1 

430 

-1.4 + 1 

36, 22 

j1544+4937 

23.226 

0.45(3) 

21(2) 

11.5(5) 


40 

5.4 

322 

-2.3 

37 

J1640+2224 

18.426 

0.21(1)* 

6.7(4)* 

15.1(9) 

450 + 200'' 

370+ 180 

2+ 1 

1472 

-2.18 + 0.12 

4, 38 

J1709+2313 

25.347 

0.80(7) 

17(2) 

3.4(6) 


30 + 5 

2.52 + 0.07 

430 

-2.1 +0.1 

35 

J1713+0747 

15.992 

0.44(4) 

9.5(8) 

11(1) 

250+ 100 

200 + 20 

36 

430 

-1.5 +0.1 

12, 4 

J1730-2304 

9.617 

0.78(6) 

9.6(8) 

130(20) 

310+ 120 

350+ 120 

43 + 6 

436 

-1.8 +0.3 

13, 4 


Notes, ^mean etTors in parentheses are only nominal uncertainties. The realistic uncertainties are estimated to be within 50% of the actual flux density values (see Sect. 3.2.2). Flux densities by 
Kuzmin & Losovsky (2001, KLOl) are at 102 or 111 MHz. Observations were done at 111 MHz. Observations were done both at 102 and 111 MHz, and it is not clear for which frequency the 
flux densities were reported. **' Profile does not show an evident scattering tail. ® Subject to substantial intrinsic variability, as pointed out by both Archibald et al. (2009) and Deller et al. (2012). 

References. (1) Lommen et al. (2000); (2) Toscano et al. (1998); (3) Navarro et al. (1995); (4) Kramer et al. (1998); (5) Burgay et al. (2006); (6) Lundgren et al. (1995); (7) Nicastro et al. (1995); 
(8) Sayer et al. (1997); (9) Archibald et al. (2009); (10) Bailes et al. (1997); (11) Wolszczan & Frail (1992); (12) Camilo (1995); (13) Lorimer et al. (1995); (14) Hessels et al. (2011); (15) Stairs et al. 
(2005); (16) Foster et al. (1991); (17) Lorimer et al. (1996); (18) Janssen et al. (2010); (19) Champion et al. (2005); (20) BoriakofF et al. (1984); (21) Fruchter et al. (1990); (22) Bates et al. (2013); 
(23) Stappers et al. (1996); (24) Cognard et al. (2011); (25) Camilo et al. (1996b); (26) Nice et al. (1993); (27) Ransom et al. (2011); (28) Lynch et al. (2013); (29) Jacoby (2005); (30) Ransom et al. 
(2014); (31) Deller et al. (2012); (32) Stovall et al. (2014a); (33) Hobbs et al. (2004a); (34) Lorimer et al. (2005); (35) Lewandowski et al. (2004); (36) Lorimer et al. (2007); (37) Bhattacharyya 
et al. (2013); (38) Kuniyoshi et al. (2015). 
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Table 4. continued. 


PSR 

DM 

IVeff 

6 

mean 

KLOl flux density 

Predicted S mean 

5i-ef 

Vi-ef 

Spectral 

Reference 


(pc cm“^) 

(ms) 

(%) 

(mJy) 

(mJy @ 102 MHz) 

(mJy) 

(mJy) 

(MHz) 

index 


J1738+0333 

33.778 

0.58(5) 

9.9(8) 

17(3) 


50 + 100 

2 

1400 

-1.4 +1 

29, 22 

J1744-1134 

3.139 

0.23(3)* 

5.6(8)* 

38(5) 

220+ 100'' 

160 + 20 

18 + 2 

436 

-1.85 + 0.08 

2,38 

J1810+1744 

39.659 

0.566(3) 

34.1(2) 

563(1) 


240 + 30 

20 

350 

-2.57 + 0.16 

14,38 

J1816+4510 

38.887 

0.69(2) 

21.7(8) 

83.5(5) 


10 

1.5 

350 

-2 

32 

J1853+1303 

30.570 

0.75(6) 

18(2) 

18(2) 


10 + 20 

0.4 + 0.2 

1400 

-1.4 +1 

15,22 

B1855+09 

13.300 

0.85(4) 

15.9(8) 

32(3) 

450 + 250 

140 + 30 

32 + 6 

430 

-1.3 +0.2 

16 

J1905+0400 

25.692 

0.7(1) 

18(3) 

21(4) 


1+2 

0.05 + 0.01 

1400 

-1.4 +1 

33,22 

J1911-1114 

30.975 

0.64(6) 

18(2) 

38(4) 

260+ 130 

470+ 110 

31+9 

400 

-2.45 + 0.06 

17,38 

J1918-0642 

26.554 

1.9(1) 

24(2) 

38(3) 


30 + 5 

5.9 + 0.6 

350 

-1.67 + 0.06 

18 

J1923+2515 

18.858 

0.37(1)* 

9.7(4)* 

17(1) 


15 + 5 

0.6 + 0.2 

820 

-1.7 

28 

B1937+21 

71.040 

0.697(6) 

44.7(4) 

370(2) 


4950 + 490 

232 + 25 

430 

-2.59 + 0.04 

16,38 

j1944+0907 

24.34 

1.48(2) 

28.6(4) 

81(2) 


20 + 20 

3.9 + 0.3 

430 

-1.4 +1 

19, 22 

B1953+29 

104.501 

1.98(5) 

32.3(8) 

38(1) 


200 + 80 

15 + 5 

430 

-2.2 +0.3 

20, 4 

B1957+20 

29.117 

0.576(6) 

35.8(4) 

155(1) 


400-1300 

^^-13 

430 

-2.79 + 0.11 

21,38 

J2019+2425 

17.203 

0.51(6) 

13(2) 

6(1) 

190+ 100 

15 + 15 

2.7 + 0.7 

430 

-1.4 +1 

26, 22 

J2043+1711 

20.710 

0.415(9)* 

17.4(4)* 

34(1) 


3 + 2 

0.8 

350 

-1.4 +1 

14, 22 

J2051-0827 

20.745 

0.48(2)* 

10.7(4)* 

43(3) 

250+ 100' 

180 + 60 

22 

436 

-1.8 +0.3 

23, 4 

J2145-0750 

8.998 

1.17(2)* 

7.3(1)* 

162(2) 

480+ 120 

470+ 130 

100 + 30 

436 

-1.33 + 0.18 

2,38 

J2214+3000 

22.557 

0.47(5) 

15(2) 

5.4(7) 


30 + 40 

2.1 

820 

-1.4 +1 

27, 22 

J2215+5135 

69.2 

0.76(1) 

29.0(4) 

53(2) 


120 + 20 

5 

350 

o OO+0.15 
*^•^^-0.23 

14,38 

J2235+1506 

18.09 

2.5(2)* 

4.2(4)* 

6(1) 

170 + 80* 

15 + 15 

3 

430 

-1.4 +1 

12, 22 

J2302+4442 

13.762 

0.64(4) 

12.4(8) 

8.8(7) 


30 + 60 

1.2+ 0.4 

1400 

-1.4 +1 

24, 22 

J2317+1439 

21.907 

0.27(1)* 

7.7(4)* 

37.9(8) 

90 + 50 

90 + 30 

19 + 5 

430 

1 

1+ 

p 

25, 4 

J2322+2057 

13.372 

0.3(2) 

7(3) 

2(1) 

75 + 40 

15 + 15 

2.8+ 1.2 

430 

-1.4 +1 

25,22 


Notes. 5mean errors in parentheses are only nominal uncertainties. The realistic uncertainties are estimated to be within 50% of the actual flux density values (see Sect. 3.2.2). Flux densities by 
Kuzmin & Losovsky (2001, KLOl) are at 102 or 111 MFIz. Observations were done at 111 MHz. Observations were done both at 102 and 111 MHz, and it is not clear for which frequency the 
flux densities were reported. Profile does not show an evident scattering tail. ® Subject to substantial intrinsic variability, as pointed out by both Archibald et al. (2009) and Deller et al. (2012). 

References. (1) Lommen et al. (2000); (2) Toscano et al. (1998); (3) Navarro et al. (1995); (4) Kramer et al. (1998); (5) Burgay et al. (2006); (6) Lundgren et al. (1995); (7) Nicastro et al. (1995); 
(8) Sayer et al. (1997); (9) Archibald et al. (2009); (10) Bailes et al. (1997); (11) Wolszczan & Frail (1992); (12) Camilo (1995); (13) Lorimer et al. (1995); (14) Hessels et al. (2011); (15) Stairs et al. 
(2005); (16) Foster et al. (1991); (17) Lorimer et al. (1996); (18) Janssen et al. (2010); (19) Champion et al. (2005); (20) Boriakoff et al. (1984); (21) Fruchter et al. (1990); (22) Bates et al. (2013); 
(23) Stappers et al. (1996); (24) Cognard et al. (2011); (25) Camilo et al. (1996b); (26) Nice et al. (1993); (27) Ransom et al. (2011); (28) Lynch et al. (2013); (29) Jacoby (2005); (30) Ransom et al. 
(2014); (31) Deller et al. (2012); (32) Stovall et al. (2014a); (33) Hobbs et al. (2004a); (34) Lorimer et al. (2005); (35) Lewandowski et al. (2004); (36) Lorimer et al. (2007); (37) Bhattacharyya 
et al. (2013); (38) Kuniyoshi et al. (2015). 
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Fig. 3. Measured mean flux densities at 110-188 MHz versus predicted 
flux densities using high-frequency data for 48 MSPs. The dotted line 
shows the loci of measured mean fluxes equal to predicted fluxes, and 
dashed lines mark the region where measured mean flux density devi¬ 
ates from predicted flux density by a factor of two. Open symbols show 
the MSPs for which spectral indices are not known, and the average 
index from Bates et al. (2013) was used, while solid symbols are for 
those with known spectral indices (see Table 4). Blue squares and green 
circles indicate MSPs with measured duty cycles > 20% and < 20%, 
respectively. Pulsars are labelled by their right ascension. 

dieted flux density comes from our poor knowledge of the spec¬ 
tral indices. However, for about another third of the MSPs, 
the measured flux densities are still lower than the predicted 
ones even including the uncertainty on their spectral indices. 
This can be explained by contributions from other factors given 
in Sect. 3.2.1. For instance, the profiles of PSRs B1937-I-21, 
B1957-I-20, J0218-I-4232, B1953-H29 are dominated by scatter¬ 
ing, so our flux density measurements are underestimated for 
them. PSR J1023-I-0038 has a very large predicted flux of up 
to ~ 1.7Jy. This could be due to significant intrinsic variability 
(Archibald et al. 2009). For others, we note that their profiles 
are quite narrow and show only a small scattering tail or none at 
all. It is thus possible that the spectra of these pulsars turnover 
within or above our observing band. For example, Kuzmin & 
Losovsky (2001) report a turnover for PSR J1012-I-5307 near 
100 MHz. We can not rule out this possibility, but more accurate 
measurements are needed at both HBA and LBA bands. Also, 



Predicted mean flux density (mjy) 


Fig. 4. Same as the Fig. 3, but predicted mean flux densities are those 
from Kuzmin & Losovsky (2001) scaled to the frequency range of 110- 
188 MHz using the spectral indices from Table 4. 


Kuniyoshi et al. (2015) recently reported that PSRs J1640-1-2224 
and J1911-1114 likely have a spectral break near 100MHz. 

Refractive scintillation could also affect our measurements, 
although not by more than a factor of ~ 1.5 (see Sect. 3.2.1). 
Another possible explanation for the discrepancy between mea¬ 
sured and predicted flux density values, is that published flux 
densities at high frequencies were over- or underestimated, most 
likely due to scintillation as was recently noted by Levin et al. 
(2013). To exclude the influence of overestimated flux densities 
at higher frequencies, we used the flux density measurements 
at frequencies 102/111 MHz from Kuzmin & Losovsky (2001) 
to predict flux densities at 150 MHz. The corresponding plot is 
shown in Fig. 4. If there is a real turnover in the spectra we 
would expect our flux density measurements to be higher than 
the predicted flux density and, hence, located in the top-left part 
of the plot. Our flux density value for PSR J1012-I-5307 is in¬ 
deed large which speaks in favour of spectral turnover for this 
MSP. The same could also apply for PSR J0034-0534, although 
its flux density is consistent with what is expected from high- 
frequency measurements, and the flux density value reported by 
Kuzmin & Losovsky (2001) could be underestimated. For PSR 
J1640-I-2224 Kuniyoshi et al. (2015) report a possible turnover 
below 100 MHz, however our flux measurement is still lower 
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than that predicted from the Kuzmin & Losovsky (2001) mea¬ 
surements at 102/111 MHz. Assuming all flux measurements are 
not affected by other factors (see Sect. 3.2.2), turnover must oc¬ 
cur sharply below 100 MHz. Similarly to the predictions from 
high-frequency data, our measured flux densities lie below the 
predicted ones for about ten out of 21 pulsars with available 
flux densities from Kuzmin & Losovsky (2001). Not ruling out 
other explanations including our own flux densities being some¬ 
what underestimated due to scattering, RISS, beam jitter or con¬ 
tribution to the system temperature from the Galactic plane in 
the FoV, the 100-MHz flux density measurements from Kuzmin 
& Losovsky (2001) could be subject to a significant bias due 
to diffractive scintillation. The decorrelation bandwidth, Avd, at 
100 MHz is of the orderof < 80 kHz (Cordes et al. 1985), and the 
receiver bandwidth, B, of Kuzmin & Losovsky (2001) was only 
160 kHz, comparable to the frequency scale of diffractive scintil¬ 
lation. This would modulate the pulsar flux density by a factor of 
VAvd/B (Popov & Soglasnov 1984). If a pulsar is weak enough 
to be only detected during the favourable ISM conditions, then 
the measured flux density will be overestimated. 

LOFAR’s observing bandwidth is much larger than the char¬ 
acteristic decorrelation bandwidth at 150 MHz (Avd < 0.2 MHz), 
so diffractive scintillation cannot affect our measurements. How¬ 
ever, to make firm flux density estimates one needs to perform a 
long-term monitoring campaign to account for RISS. 

3.3. LB A detections 

We also observed nine bright MSPs with the LOFAR LBAs 
in the frequency range 15-93 MHz. We detected three MSPs, 
namely PSRs J0030H-0451, J0034-0534, and J2145-0750. All 
detected and non-detected MSPs in LBA observations are listed 
in the last column of Table 2. Only PSR J2317-I-1439 was ob¬ 
served for 20 min, for the others the integration time was 1 h. 

It is currently unclear which effect plays the dominant role 
in the non-detections of MSPs below 100 MHz: scattering, in¬ 
creased sky temperature, or smaller flux densities due to a 
possible turnover in their spectra. However, using the empiri¬ 
cal relationship from Bhat et al. (2004) for DM = lOpccm^^ 
one would expect the scattering time to be about 0.37 ms at 
57 MHz, or about ~ 0.07 cycles or less for PSRs J1012H-5307, 
J1022H-1001, J1024-0719, and B1257-H12. Hence, scattering is 
most likely not the primary reason for the non-detections of 
these four MSPs. From HBA profiles of these pulsars we deter¬ 
mined rough upper limits on the Tjcat at 57 MHz to be 9, 6.4, 28, 
and 9 ms, respectively. The limits do not contradict the predic¬ 
tion of Tscat from the Bhat et al. (2004) relationship. However, 
although rough, our limits still give the possibility for scatter¬ 
ing to play a role in the non-detections of PSRs J1012H-5307, 
J1024-0719, and B1257-1-12, as the uncertainty in the Bhat et al. 
(2004) relationship can be up to two orders of magnitude. For 
PSR J1022-I-1001 the spin period is about three times longer than 
our rough limit on Tscat, thus most likely other factors contribute 
to the non-detection of this MSP - more so than scattering. 

These four non-detected MSPs are located at high Galactic 
latitudes similar to the three MSPs detected with the LBAs, and 
the sky temperatures in these directions are very similar. How¬ 
ever, the spectral indices for two of the detected MSPs, PSRs 
J0030H-0451 and J0034-0534 (-2.2 + 0.2 and -2.6 + 1, respec¬ 
tively, see Table 4; approximately equivalent to the spectral in¬ 
dex of Galactic synchrotron radiation), are steeper than for these 
four non-detected MSPs. Hence, the influence of the increased 
background temperature at low radio frequencies for PSRs 
J0030H-0451 and J0034-0534 is much smaller than for the non- 


detected MSPs with shallower spectral indices. Alternatively, 
these four non-detected MSPs might have a turnover in their 
radio spectra at frequencies < 120MHz (for PSR J1012-I-5307 
Kuzmin & Losovsky (2001) reported a possible turnover in its 
spectrum near 100 MHz). 

For PSR J2317H-1439 the predicted scattering time at 
57 MHz is about 3 ms, comparable to the pulse period. The rough 
upper limit on Tjcat from the HBA profile is somewhat lower, 
about 2.5 ms. Thus, similarly to the other four non-detected 
MSPs, scattering is not the dominant factor in its non-detection. 

For the black widow PSR J1810-1-1744, with a DM of 
39.7pccm“^, the scattering time is already four times longer 
than the spin period, and the profile is scattered out (Archibald 
et al. 2014). One can already see profile evolution within the 
HBA band for this MSP in Fig. 9 (left), where the profile at the 
low-frequency edge is notably more scattered than at the high- 
frequency end of the HBA band. 

For completeness, we should also consider the influence of 
iiorb/c smearing on LBA profiles of binary MSPs as was dis¬ 
cussed in Sect. 3.1.1, as this effect becomes even more sig¬ 
nificant at lower observing frequencies. PSRs J0030-t-0451 and 
J1024-0719 are isolated, so this effect is not relevant for them. 
For the remaining seven MSPs we calculated the extra smear¬ 
ing of the profile due to this effect in a 3.9-MHz frequency 
band near a frequency of 58 MHz, where we have our high¬ 
est sensitivity. For four out of seven MSPs (PSRs J1012-I-5307, 
J1022-H1001, B1257-bl2, and J2145-0750) this effect is minor 
with less than 0.01 cycles of extra broadening of the profiles. 
For PSR J2317-1-1439, which we did not detect, the broadening is 
0.04 cycles, still quite small, so we should have been able to de¬ 
tect it regardless of this effect. In contrast, for PSR J1810-I-1744 
the effect is quite significant producing additional profile scatter¬ 
ing of about 0.1 cycles. However, as we already showed, the ex¬ 
pected broadening due to the scattering itself is 100 times larger. 

Figure 5 shows the profile evolution for the three detected 
MSPs. The profile on top of each panel is the average profile 
from 38-77 MHz. Among these three, PSR J0034-0534 has the 
largest DM of 13.8pccm“^ and shows significant profile evo¬ 
lution with its profile clearly being scattered towards lower fre¬ 
quencies. The profile shape in the total LBA band is largely de¬ 
termined by scattering at the lower part of the band, although the 
profile in the higher-frequency sub-band at 73.3 MHz resembles 
the profile in the HBA band (see Fig. 1). At this observing fre¬ 
quency we measured the scattering time to be about 440 p.s by 
fitting the trailing part of the profile with an exponential function. 
The extra profile broadening due to the Corb/c effect is about 0.05 
cycles for a 7.8-MHz sub-band around 57.7 MHz. It does scatter 
the profile but not at a large enough level to scatter it out com¬ 
pletely for the lowest sub-band shown in Fig. 5 (middle). There¬ 
fore, the profile broadening must be dominated by scattering in 
the ISM. 

The profiles of the other two MSPs, PSRs J0030-I-0451 
and J2145-0750, do not show significant evolution in com¬ 
parison to PSR J0034-0534. They have smaller DMs (4.3, 
9.0, 13.8pccm-3 for PSRs J0030-H0451, J2145-0750, and 
J0034-0534, respectively), and the scattering conditions must 
be temperate towards them. We measured the scattering time 
to be about 300 p.s for PSR J0030-I-0451 at 42.1 MHz and less 
than 700 qs for PSR J2145-0750 at 57.7 MHz. Although there 
is no significant broadening of the profile for PSR J2145-0750 
its profile almost vanishes at 42 MHz. For PSR J2145-0750 the 
sky temperature is four times higher at 42 MHz than at 73 MHz. 
The LOFAR LBA bandpass response is not uniform and peaks 
at approximately 58 MHz (see van Haarlem et al. 2013, Fig. 19). 
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Fig. 5. Profile evolution of the three MSPs detected in the LBA frequency range of 38.2-77.2 MHz. The grey profile on top of each panel is the 
average profile within the 39-MHz wide band. All profiles are from 1-h LBA observations. There are 64 bins in the profiles of PSRs J0030+0451 
and J0034-0534 and 128 bins in the profiles of PSR J2145-0750. The numbers on the left in each panel show the central frequency in MHz of 
the corresponding 7.8-MHz sub-band. 


However, we do not think that lower sensitivity is the reason 
for the much weaker profile at 42 MHz. For PSR J0030-P0451 
the change in the sky temperature and bandpass response is 
the same, but it has a clearly detected profile at 42 MHz, only 
slightly weaker than the profile at 58 MHz. Hence, we believe 
that this weakening is due to a turnover in the spectrum of 
PSR J2145-0750. Our assumption qualitatively agrees with the 
recent conclusion by Dowell et al. (2013), who found that the 
radio spectrum of PSR J2145-0750 is best fit with a model with 
spectral curvature and a roll-over frequency of 730 MHz using 
data from the LWA. Recently, Kuniyoshi et al. (2015) also re¬ 
ported on a spectral turnover between 100 and 400 MHz based 
on archival and LWA data. 

Stovall et al. (2014b) also report on the detection of the same 
three MSPs with the LWA between 30-88 MHz from their 6-h 
observations of PSR J0030-b0451 and 8-h observations of PSRs 
J0034-0534 and J2145-0750. 

4. Discussion 

4.1. Detectability of MSPs 

In order to explore the degree to which propagation in the ISM 
affects the detectability of the MSPs presented here, we plotted 
the ratio of the predicted Tscat at 150 MHz, over the pulsar pe¬ 
riod, P, on DM and S’mean at 110-188 MHz. This dependence is 
shown in Fig. 6 for both detected and non-detected MSPs, with 
the green colour indicating the 5mean- For non-detected MSPs, 
the green intensity shows the upper limit on 5 mean as given in 


Table 3 (Col. 10). For Tscat we used the NE2001 Galactic free 
electron density model (Cordes & Lazio 2002) and assumed a 
Kolmogorov spectrum of electron density inhomogeneities (see 
Table 5, last column). 

The histogram on top of Fig. 6 presents the number of de¬ 
tected and non-detected MSPs as a function of DM. We de¬ 
tected 25 out of 28 observed MSPs with DMs < 20pccm^^. 
On the high-DM side, we did not detect two MSPs with DM 
> 110 pc cm“^ that we tried, namely PSR J1949-1-3106 with DM 
= 164pccm-3 and PSR J1903-1-0327 with DM = 297 pc cm ^ 
On the other hand, one of the detected MSPs, PSR B1953-b29, 
has a DM of 104.5 pc cm“^. For the DM range of 20-100pc cm“^ 
we have almost parity between detected (22) and non-detected 
MSPs (23). 

There is no clear dependence of the MSP detectability on the 
predicted scattering measure. If it were the case, we would see 
the majority of non-detected MSPs clustering close to or above 
the line Tscat/P = 1, where a pulse profile would be completely 
washed out by the scattering. However, most of our non-detected 
MSPs lie below Tscat/P = 1 line populating the DM range of 
~20-50pccm“^. The dependence of Tjcat/P on DM follows the 
empirical relationship between DM and Tscat provided by Bhat 
et al. (2004). However, there is more than two orders of magni¬ 
tude of scatter in this relationship, hence its predictive power for 
individual pulsars is very limited. Therefore, it is certainly pos¬ 
sible that for some of the non-detected MSPs the scattering con¬ 
ditions are much worse than we can infer from their DM (and/or 
their detectability depends on other factors). 
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Fig. 6. MSP detectability, or dependence of the ratio of predicted scattering time at 150 MHz over the pulsar period, Tjcat/P. on dispersion measure, 
DM. The green colour map corresponds to the measured mean flux density, Smean, at 110-188 MHz. Non-detected MSPs are shown with brown 
squares and their green intensity level corresponds to the measured upper limit of the mean flux density at 110-188 MHz. Scattering times at 
150 MHz were calculated from the NE2001 Galactic free electron density model (Cordes & Lazio 2002) assuming a Kolmogorov spectrum of 
electron density inhomogeneities. The dashed line shows the empirical dependence between Tjcat and DM at 150 MHz from Bhat et al. (2004), 
where Tjcat is arbitrarily normalised by P = 4.075 ms, which is the median value of MSP periods in our sample of detected and non-detected 
MSPs. The top subplot shows the histograms of DMs for detected and non-detected MSPs. The inset plot is the same as the main plot but only 
for the non-detected PSR J1903-1-0327, which did not fit into the range of the main plot. PSR J1400-1438 was recently detected with LOFAR 
using the updated coordinates (J. Swiggum, private communication) derived from over three years of timing with the GBT The previously used 
catalogued position was off by about 6', which is about two times larger than the full-width at half-maximum of LOFAR’s full-core tied-array 
beam at 150 MHz. 


Another factor that can play an equal or even larger role 
than scattering in the MSP detectability is our lack of sen¬ 
sitivity resulting from the increased sky temperature back¬ 
ground at 150 MHz. In Fig. 6, all five detected MSPs with 
DM > 40 pc cm ^ have 5 mean above 27 mJy, larger than the 
average value for the rest of the detected MSPs. Also, except 
for the double pulsar, PSR J0737-3039A, four others pulsars 

- PSRs J0218-H4232, B1937-H21, B1953H-29, and J2215-H5135 

- have significantly scattered profiles, and their 5 mean most 
likely is underestimated. For PSRs J0218-I-4232, B1937-H21, and 
B1953-1-29 spectral indices derived from high-frequency data 
are -2.8, -2.3, and -2.2, respectively, significantly larger than 
the average for the pulsar population and comparable with the 
spectral index of the synchrotron background radiation. In fact, 
PSRs B1937-b21 and J0218 -1-4232 were first recognised as steep- 
spectrum, highly linearly polarised, compact radio sources in in¬ 


terferometric imaging data (Erickson 1980; Navarro et al. 1995), 
which is what led to their discovery. Most of the other detected 
MSPs have shallower spectral indices. Together with the pos¬ 
sible turnover in their spectra and scattering, this could hinder 
MSP detection at 150 MHz. For instance, PSR J1713-1-0747 is 
very bright (S/N ~ 4000) at 1.4 GHz (see e.g. Dolch et al. 2014). 
However, the LOFAR detection was very weak, with 5 mean of 
only about 11 mJy (S/N ~ 8), and the profile is not very scat¬ 
tered (Tscat ~ 120 ps). Thus, in this case the main reason for a 
weak detection is not due to the scattering, but rather the spec¬ 
tral index of roughly -1.5 being smaller than for the synchrotron 
background radiation, or the possible turnover in its radio spec¬ 
trum between 200 and 1400 MHz, or both. 

Three of our non-detected MSPs, PSRs J0613-0200, 
B1620-26, and J2229-I-2643, were previously reported to be 
detected at 102/111 MHz by Kuzmin & Losovsky (2001). As 
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Fig. 7. Dispersion measure versus the MSP duty cycle, d. Green circles show MSPs with profiles which do not show an evident scattering tail - 
see S values marked by an asterisk (*) in Table 4. Light yellow squares show all other MSPs. Errors in DMs including dDMo^b (see Sect. 4.3) are 
smaller than the size of the marker for all MSPs. For MSPs without error bars, the uncertainty on d is also less than the size of the marker. The top 
subplot shows the histogram of duty cycles for all detected MSPs. 


we discussed in Sect. 3.2.3, their detection and flux measure¬ 
ments could be significantly biased due to diffractive scintilla¬ 
tion. PSR J0613-0200 was observed by Stappers et al. (2008) 
with the WSRT at 150 MHz and they did not detect it. Also, for 
these and other non-detected MSPs additional factors discussed 
in Sect. 3.2.1 could affect their detection, most notably the iono¬ 
spheric beam jitter and RISS. 


4.2. Profile widths 

Table 4 lists measured values of Wefr and S in Cols. 3 and 4 for the 
total observing bandwidth of 110-188MHz. At low observing 
frequencies scattering broadens the intrinsic pulse profiles, and 
our measured VTeff is essentially affected by scattering for many 
of the observed MSP profiles. For MSPs with narrow, visibly 
unscattered profiles, Wefr and 6 in Table 4 are marked with an 
asterisk (*). Figure 7 shows 6 vs DM for detected MSPs. Almost 
all MSPs with 6 > 20% show profiles widened by scattering 
(yellow circles). Also, between them there is a trend of d being 
larger with an increasing DM. The only exception here is PSR 
J0034-0534 which has a very broad profile (S ^ 29%), which 


mostly likely is genuinely broad by comparing with the WSRT 
profile at 376 MHz (see Fig. 10). 

For the majority of the detected MSPs, d < 20%. Among 
them, MSPs with low DMs < 20pccm“^ have sharp narrow 
profiles and are weakly (if at all) affected by scattering (green 
circles). Others, at moderately high DMs, have profiles that are 
strongly affected by scattering and/or intrinsically weak at 110- 
188 MHz. 

Table 4 gives nominal 1-cr errors on Wefj and b values. The 
actual error for Wefj can be up to two times larger for the majority 
of binary MSPs depending on the ratio Uorb/c (see Sect. 3.1.1). 
The actual error on 6 could be larger by 0.1 %, and does not affect 
our conclusions on MSP widths. 

4.3. DM variations 

It is known that, on time scales of months to years, DM and Tjcat 
towards some pulsars change due to the varying line-of-sight 
through the ISM (e.g. Hobbs et al. 2004b). For instance, Hem- 
berger & Stinebring (2008) showed that detectable variations of 
the scattering measure towards PSR B1737H-13 occurred roughly 
on a time scale of a fortnight for observations in the 1.1-1.5-GHz 
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range. Keith et al. (2013) reported on long-term DM variations 
for MSPs on time scales of a few months. Correction for such 
DM variations is crucial for high-precision timing. 

To measure DM at the epochs of the individual observations 
we split the observing band into four-five sub-bands and for 
each of them generated a single TOA for the whole observation 
by cross-correlation with a profile template. The template was 
based on the average pulse profile of all available observations 
in the total HBA band. We then used the pulsar timing package 
Tempo2^^ (Hobbs et al. 2006) to fit the generated TOAs only for 
DM. In Table 5 we present the weighted average, DMhba, of 
our DM measurements (Col. 5) for each of the individual ob¬ 
serving epochs during the given time span (Col. 7). The refer¬ 
ence epoch that corresponds to the middle of the observing time 
span is given in Col. 6. The number of observations, Aobs, per¬ 
formed within a given span is listed in Col. 8. For 11 out of 48 
detected MSPs we carried out only one exploratory observation, 
therefore the presented DMs are individual DM measurements 
for the epoch of the observation. We calculated the DM offset, 
ADM = DMref - DMhba, between the catalogue DM value, 
DMref, and LOFAR’s DM measurement, DMhba (Col. 9), with 
the uncertainty, cr^DM, being the quadrature sum of uncertainties 
of both DMref and DMhba- 

Column 11 of Table 5 lists measured DMs for the best in¬ 
dividual observations. The uncertainties for most of the MSPs 
are lower than the errors on DMref, owing to the low observ¬ 
ing frequency and large fractional bandwidth of LOFAR. The 
uncertainty on the weighted average DM value, DMhba, on the 
other hand is larger than the uncertainty of the individual mea¬ 
surements for the majority of the MSPs. This is most likely ev¬ 
idence for short-term DM variations. However, care should be 
taken to account for artificial orbital DM variations for binary 
MSPs as was discussed in Sect. 3.1.1. These DM variations, 
bDMorb = forb/c X DM, are larger for higher-DM pulsars in com¬ 
pact binary orbits. We list bDMorb for all binary MSPs in Col. 10 
of Table 5. For Uorb we used the maximum orbital velocity of the 
pulsar at periastron, lip/c = (27rA/Pb)[(l+e)/(l-e)]'^^ wherexis 
the projected semi-major axis, e is the eccentricity, and Pb is the 
orbital period of the pulsar. In our binary MSP sample IdDMoitl 
is highest for the double pulsar J0737-3039A (0.054pccm“^), 
and lowest for PSR B1257-I-12 (< lO^^pccm”^). The influence 
of bDMorb on the weighted average DMhba and ADM should be 
even smaller, as the actual dDMorb depends on Vorb at the epoch 
of observation. 

Armstrong et al. (1995) have shown that the structure of the 
ISM has a steep red power spectrum of electron density that 
results in DM variations becoming more prominent at longer 
time scales. Hence, one would expect to see larger DM offsets 
for larger epoch offsets. In Fig. 8 we plotted DM offset, ADM, 
against the epoch offset between the reference catalogue epoch 
and LOFAR epoch of DM measurements. To account for differ¬ 
ent uncertainties of DM measurements for different MSPs and, 
more importantly, to exclude the influence from dDMorb for bi¬ 
nary MSPs, we normalised ADM by the quadrature sum of its 
uncertainty, cr adm, and bDMorb- Note, that for the average DM 
value over a given time span bDMorb should be also taken as the 
average contribution from all (Nobs) observations. We will, how¬ 
ever, use the most conservative estimate on dDMorb to avoid any 
influence on our conclusions. 

For about 25% of the MSPs, ADM is more than three times 
its DM uncertainty including dDMoit as well. The rest of the DM 
offsets lie within +3crADM area including all MSPs for which the 

https://bitbucket.org/mkeith/tempo2 


bDMoib is larger than cr adm (brown error bars). We do not see 
any obvious trends for ADM between binary and isolated MSPs. 
Finally, we also do not see any trends for ADM to become larger 
for larger epoch offsets. There are six MSPs with large ADM 
within five years from the reference DM epoch. This discrepancy 
in DM may be due to large secular DM drift in the correspond¬ 
ing line of sights. Another possibility for large DM differences 
at such short-term scales may reflect the fact that at low frequen¬ 
cies the radio emission probes somewhat different regions of in¬ 
terstellar space, thus providing a slightly different DM (Cordes 
et al. 2015). Multifrequency simultaneous observations and DM 
measurements are needed to confirm this point conclusively or 
rule it out. 

Keith et al. (2013) presented DM variations for 20 MSPs 
from the Parkes Pulsar Timing Array (PPTA); eight sources 
overlap with our sample. Unfortunately, observing time spans do 
not overlap, with the last measurement in their paper being for 
MJD 55800. Hence, we cannot make a direct comparison of DM 
offsets. However, a few of the overlapping MSPs show large reg¬ 
ular DM drifts, in particular PSRs J1730-2304 and B1937-F21, 
so we can try to extrapolate the DM offsets at our observing 
epochs for at least these two MSPs. For PSR J1730-2304 the 
extrapolation to the reference epoch of MJD 56474 gives ADM 
of about 4.5x10“^ pc cm“^, comparable with our measured offset 
of 0.006 pc cm“^. We see similar agreement for PSRB1937H-21, 
where DM shows a secular decrease with an expected value 
at the reference epoch of MJD 56385 of ~ 5 x lO^^pccm^^ 
smaller than the DMref. This is somewhat smaller than our off¬ 
set of 0.016 pc cm“^, but the amount of DM drift in Keith et al. 
(2013) for this MSP does show significant variations (see also 
Kaspi et al. 1994). One of the three LOFAR observations of PSR 
B1937-I-21 at MJD 56258.75 also overlaps with the intense mul¬ 
tifrequency observing campaign of this pulsar presented in Yi 
et al. (2014). Our DM measurement for this observing epoch is 
about 5 X 10“^ pc cm“^ larger than theirs. In their timing analysis 
they accounted for the broadening of the profile due to scattering 
in each frequency channel across the 310-390MHz band, which 
likely explains the observed DM offset between our measure¬ 
ments. 

Precise LOFAR DM measurements on time scales of two- 
four weeks could provide valuable data for DM corrections of 
higher-frequency pulsar timing data, especially for the MSPs 
that do not show variations within the uncertainties of their 
measurements at higher frequencies, e.g. PSRs J1022-I-1001, 
J1713-I-0747, and J2145-0750. We have already started regu¬ 
lar monthly observations of MSPs with LOFAR and we also 
perform observations with shorter cadence using international 
LOFAR stations. With these regular timing observations we will 
measure DM values and provide them to the broader pulsar com¬ 
munity for use in ISM studies and high-precision timing. A de¬ 
tailed analysis of DM variations will be presented elsewhere. 

4.4. Profile variations 

In the previous Section we assumed there is no profile varia¬ 
tion within the observing band. However, at LOFAR frequencies 
this is usually not the case as there are several factors that can 
affect the observed profile shape within the band, such as the 
intrinsic profile evolution and profile smearing due to scatter¬ 
ing. For pulsar timing, normally a single profile template is used 
for cross-correlation with the observed profile to get a TOA. If 
the profile at the low-frequency edge of the band is significantly 
scattered, it will artificially increase the uncertainty at the lower 
end of the band and produce a later TOA, which can be wrongly 
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Table 5. Dispersion measures of detected MSPs in the frequency range 110-188 MHz. 


PSR 

Catalogue 

DMref Epoch 

(pccm-3) (MJD) 

Ref 

DMhba 
( pc cm“^) 

LOFAR 

Epoch 

(MJD) 

Span 

(days) 

A^obs 

ADM 
(pc cm“^) 

ADMo,-b 
(10-3 pc 
cm-3) 

DMxempo2 
(pc cm“^) 

DMpditip 

(pc cm“^) 

^scat 

(ps) 

J0030+0451 

4.333(1) 

50984.40 

1 

4.3326(2) 

56371.00 

154.6 

6 

+0.0004(10) 

— 

4.332613(14) 

4.3326(1) 

1.7 

J0034-0534 

13.76517(4) 

50690.00 

2 

13.76508(17) 

56367.52 

161.6 

6 

+0.00009(17) 

0.906 

13.765034(7) 

13.76498(8) 

60 

J0214+5222 

22.0354(34) 

55974.00 

6 

22.03644(24) 

56716.60 

139.6 

6 

-0.0010(34) 

0.549 

22.0361(1) 

22.035(2) 

310 

J0218+4232 

61.252(5) 

50864.00 

3 

61.23889(32) 

56473.30 

— 

1 

+0.013(5) 

4.357 

61.2389(3) 

61.2395(3) 

2430 

J0337+1715 

21.3162(3) 

55920.00 

4 

21.319(1)* 

56512.23 

— 

1 

-0.0028(10) 

1.159 

— 

21.319(1) 

290 

J0407+1607 

35.65(2) 

52799.00 

20 

35.61091(17) 

56790.52 

— 

1 

+0.04(2) 

0.413 

35.61091(17) 

35.6105(9) 

780 

J0621+1002 

36.6010(6) 

50944.00 

5 

36.522(8) 

56375.26 

172.5 

6 

+0.079(8) 

3.859 

36.5349(81) 

36.5377(63) 

1110 

J0636+5129 

11.10598(6) 

56307.00 

6 

11.1068(31) 

56700.88 

105.7 

2 

-0.0008(31) 

0.109 

11.10647(21) 

11.1064(3) 

80 

J0645+5158 

18.247536(9) 

56143.00 

6 

18.24847(12) 

56552.87 

461.7 

9 

-0.00093(12) 

— 

18.24846(15) 

18.2485(6) 

260 

J0737-3039A 

48.920(5) 

53156.00 

7 

48.9179(55) 

56321.94 

— 

1 

+0.002(7) 

53.73 

48.9179(55) 

48.915(12) 

0.4 

J0751+1807 

30.2489(3) 

51800.00 

8 

30.244(1) 

56288.03 

16.0 

2 

+0.005(1) 

3.316 

30.2434(5) 

30.2425(15) 

590 

J1012+5307 

9.0233(2) 

50700.00 

9 

9.02433(8) 

56357.49 

154.6 

6 

-0.00103(22) 

0.631 

9.024355(76) 

9.0245(2) 

40 

J1022+1001 

10.2521(1) 

53589.00 

10 

10.2530(24) 

56357.46 

154.6 

5 

-0.0009(24) 

1.601 

10.25327(16) 

10.2533(1) 

50 

J1023+0038 

14.325(10) 

54802.00 

11 

14.3309(2) 

56308.15 

14.0 

2 

-0.006(10) 

1.805 

14.3311(2) 

14.3315(1) 

90 

J1024-0719 

6.48520(8) 

53000.00 

10 

6.48447(14) 

56357.46 

154.6 

4 

+0.00073(16) 

— 

6.48446(25) 

6.4846(4) 

5.6 

J1038+0032 

26.59(20) 

53000.00 

30 

26.312(2) 

56782.80 

— 

1 

+0.28(20) 

— 

26.3121(22) 

26.3069(55) 

330 

J1231-1411 

8.090(1) 

55100.00 

27 

8.09100(54) 

56789.86 

— 

1 

-0.001(1) 

0.646 

8.09100(54) 

8.0908(11) 

20 

B1257+12 

10.16550(3) 

49750.00 

12 

10.15397(13) 

56366.55 

172.5 

6 

+0.0115(1) 

< 10-3 

10.15405(11) 

10.1542(1) 

20 

J1453+1902 

14.049(4) 

53337.00 

31 

14.0569(16)* 

56779.99 

— 

1 

-0.0079(43) 

— 

— 

14.0569(16) 

90 

J1544+4937 

23.2258(1) 

56007.00 

32 

23.227(1) 

56780.00 

— 

1 

-0.0012(10) 

0.460 

23.227(1) 

23.2277(5) 

240 

J1640+2224 

18.4260(8) 

51700.00 

13 

18.42765(12) 

56366.68 

172.6 

7 

-0.00165(81) 

0.423 

18.427656(5) 

18.4280(2) 

170 

J1709+2313 

25.3474(2) 

51145.00 

33 

25.3416(24)* 

56964.53 

— 

1 

+0.0058(24) 

1.241 

— 

25.3416(24) 

320 

J1713+0747 

15.9915(2) 

52000.00 

14 

15.9907(6) 

56424.60 

286.3 

10 

+0.0008(6) 

0.555 

15.99044(35) 

15.9909(5) 

120 

J1730-2304 

9.617(2) 

53300.00 

14 

9.62314(4) 

56473.94 

33.9 

2 

-0.006(2) 

— 

9.6232(4) 

9.6233(6) 

50 


Notes. Nobs is the number of observations within the given span (Col. 7) used for calculation of DMhbaI ADM is the DM offset between the reference catalogue DM value, DM^r, and LOFAR’s 
DM measurement, DMhba. i-e. ADM = DM^r - DMhbaI dDMo^i, = forb/c x DM is conservative upper limit on “artificial” DM orbital vatiations for binary MSPs (see text); DMxcmpo 2 and DMp;jn,p 
are DM measurements for the best individual observation using Tempo2 and pditip, respectively; Tscat is the predicted scattering time at 150 MHz from the NE2001 model (Cordes & Lazio 2002) 
scaled from 1 GHz using a Kolmogorov spectrum of electron density inhomogeneities. Calculated using pdmp. 

References. (1) Abdo et al. (2009); (2) Abdo et al. (2010); (3) Hobbs et al. (2004b); (4) Ransom et al. (2014); (5) Splaver et al. (2002); (6) Stovall et al. (2014a); (7) Kramer et al. (2006b); (8) Nice 
et al. (2005); (9) Lange et al. (2001); (10) Hotan et al. (2006); (11) Archibald et al. (2009); (12) Konacki & Wolszczan (2003); (13) Ldhmer et al. (2005); (14) Verbiest et al. (2009); (15) Jacoby 
(2005); (16) Hessels et al. (2011); (17) Gonzalez et al. (2011); (18) Toscano et al. (1999); (19) Janssen et al. (2010); (20) Lynch et al. (2013); (21) Cognard et al. (1995); (22) Champion et al. (2005); 
(23) Arzoumanian et al. (1994); (24) Nice et al. (1993); (25) Guillemot et al. (2012); (26) Doroshenko et al. (2001); (27) Ransom et al. (2011); (28) Camilo et al. (1993); (29) Cognard et al. (2011); 
(30) Burgay et al. (2006); (31) Lorimer et al. (2007); (32) Bhattacharyya et al. (2013); (33) Lewandowski et al. (2004); (34) J. Hessels (private communication). 
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Table 5. continued. 


PSR 

Catalogue 

DMref Epoch 

(pccm-3) (MJD) 

Ref 

DMhba 
( pc cm“^) 

LOFAR 

Epoch 

(MJD) 

Span 

(days) 

^obs 

ADM 
(pc cm“^) 

ADMorb 
(10-3 pc 
cm-3) 

DMxempo2 
(pc cm“^) 

DMpdmp 
(pc cm“^) 

^scat 

(ps) 

J1738+0333 

33.778(9) 

52500.00 

15 

33.7754(11) 

56369.21 

175.5 

7 

+0.0026(91) 

2.378 

33.77481(44) 

33.7747(13) 

640 

J1744-1134 

3.13908(4) 

53742.00 

10 

3.13808(14) 

56353.78 

144.6 

6 

+0.0010(1) 

— 

3.13815(8) 

3.1383(2) 

<1 

J1810+1744 

39.659298(52) 

55520.78 

34 

39.66041(28) 

56349.31 

135.6 

6 

-0.0011(3) 

1.828 

39.6600(2) 

39.66050(3) 

620 

J1816+4510 

38.8874(4) 

56047.00 

6 

38.88972(13) 

56713.83 

136.7 

5 

-0.0023(4) 

4.666 

38.8897(1) 

38.8905(2) 

580 

J1853+1303 

30.5701(6) 

54000.00 

17 

30.57121(17) 

56384.24 

145.6 

2 

-0.0011(6) 

0.784 

30.57140(46) 

30.5691(9) 

130 

B1855+09 

13.300(4) 

50481.00 

14 

13.30482(89) 

56366.29 

169.5 

7 

-0.0048(41) 

0.724 

13.3047(4) 

13.3042(15) 

100 

J1905+0400 

25.6923(12) 

53700.00 

17 

25.68947(24) 

56964.66 

— 

1 

+0.0028(12) 

— 

25.68947(24) 

25.6885(38) 

130 

J1911-1114 

30.9750(9) 

50458.00 

18 

30.96909(69) 

56368.80 

164.5 

5 

+0.0059(11) 

1.462 

30.96904(43) 

30.9690(6) 

560 

J1918-0642 

26.554(10) 

53400.00 

19 

26.590(6) 

56388.75 

204.4 

3 

-0.036(12) 

1.477 

26.5913(11) 

26.5952(27) 

310 

J1923+2515 

18.85766(19) 

55322.00 

20 

18.85561(5) 

56387.77 

138.6 

2 

+0.00205(20) 

— 

18.85567(12) 

18.8557(2) 

40 

B1937+21 

71.0398(2) 

47899.50 

21 

71.024(1) 

56385.37 

253.2 

3 

+0.016(1) 

— 

71.02373(55) 

71.0252(1) 

2960 

J1944+0907 

24.34(2) 

52913.00 

22 

24.36079(38) 

56384.27 

145.6 

2 

-0.02(2) 

— 

24.3611(6) 

24.3611(7) 

80 

B1953+29 

104.501(3) 

54500.00 

17 

104.5263(26) 

56384.27 

145.6 

2 

-0.025(4) 

2.035 

104.523(1) 

104.5202(25) 

8100 

B1957+20 

29.1168(7) 

48196.00 

23 

29.12151(34) 

56322.45 

71.8 

4 

-0.0047(8) 

0.494 

29.12143(13) 

29.12130(12) 

110 

J2019+2425 

17.203(5) 

50000.00 

24 

17.1959(14) 

56401.26 

179.5 

3 

+0.007(5) 

0.634 

17.19544(84) 

17.1952(7) 

40 

J2043+1711 

20.70987(3) 

55400.00 

25 

20.7121(1) 

56384.33 

145.6 

3 

-0.0022(1) 

1.650 

20.71232(19) 

20.7124(3) 

50 

J2051-0827 

20.7449(4) 

51000.00 

26 

20.7299(17) 

56387.81 

138.6 

2 

+0.015(2) 

0.686 

20.73016(17) 

20.73120(35) 

150 

J2145-0750 

8.9977(14) 

53040.00 

14 

9.00401(13) 

56367.40 

161.6 

7 

-0.0063(14) 

0.973 

9.00401(4) 

9.00380(13) 

10 

J2214+3000 

22.557(1) 

55100.00 

27 

22.5448(22) 

56404.82 

172.5 

3 

+0.0122(24) 

0.233 

22.5431(42) 

22.5507(9) 

90 

J2215+5135 

69.2 

55135(30) 

16 

69.1976(15) 

56404.84 

172.5 

2 

+0.0024(15) 

13.14 

69.1959(7) 

69.1972(2) 

3010 

J2235+1506 

18.09(5) 

49250.00 

28 

18.10352(56) 

56521.02 

— 

1 

-0.014(50) 

— 

18.10352(56) 

18.1037(20) 

90 

J2302+4442 

13.762(6) 

55000.00 

29 

13.7268(1) 

56385.90 

148.6 

2 

+0.035(6) 

0.409 

13.72662(34) 

13.7274(5) 

20 

J2317+1439 

21.907(3) 

49300.00 

28 

21.89873(14) 

56367.48 

161.5 

6 

+0.008(3) 

1.499 

21.89876(26) 

21.8991(1) 

60 

J2322+2057 

13.372(5) 

48900.00 

24 

13.38574(25) 

56475.68 

30.9 

2 

-0.014(5) 

— 

13.3853(8) 

13.3859(9) 

30 


Notes. Nobs is the number of observations within the given span (Col. 7) used for calculation of DMhba; ADM is the DM offset between the reference catalogue DM value, DMref, and LOFAR’s 
DM measurement, DMhba. i-e- ADM = DM^f - DMhbaI dDMo^i, = Vorb/c x DM is conservative upper limit on “artificial” DM orbital vatiations for binary MSPs (see text); DMxempo 2 and DMpjj^ip 
are DM measurements for the best individual observation using Tempo2 and pditip, respectively; Tscat is the predicted scattering time at 150 MHz from the NE2001 model (Cordes & Lazio 2002) 
scaled from 1 GHz using a Kolmogorov spectrum of electron density inhomogeneities. Calculated using pdmp. 

References. (1) Abdo et al. (2009); (2) Abdo et al. (2010); (3) Hobbs et al. (2004b); (4) Ransom et al. (2014); (5) Splaver et al. (2002); (6) Stovall et al. (2014a); (7) Kramer et al. (2006b); (8) Nice 
et al. (2005); (9) Lange et al. (2001); (10) Hotan et al. (2006); (11) Archibald et al. (2009); (12) Konacki & Wolszczan (2003); (13) Ldhmer et al. (2005); (14) Verbiest et al. (2009); (15) lacoby 
(2005); (16) Hessels et al. (2011); (17) Gonzalez et al. (2011); (18) Toscano et al. (1999); (19) lanssen et al. (2010); (20) Lynch et al. (2013); (21) Cognard et al. (1995); (22) Champion et al. (2005); 
(23) Arzoumanian et al. (1994); (24) Nice et al. (1993); (25) Guillemot et al. (2012); (26) Doroshenko et al. (2001); (27) Ransom et al. (2011); (28) Camilo et al. (1993); (29) Cognard et al. (2011); 
(30) Burgay et al. (2006); (31) Lorimer et al. (2007); (32) Bhattacharyya et al. (2013); (33) Lewandowski et al. (2004); (34) J. Hessels (private communication). 
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Fig. 8. Top. Histogram of normalised DM offsets, ADM = DM^ef - DMhba- Values of ADM are normalised by the quadrature sum of DM offset 
uncertainty, cr adm, and dDMorb (see text). Bottom. Dependence of normalised ADM on the epoch offset between reference catalogue epoch and 
LOFAR epoch of DM measurements. Isolated MSPs are shown with open circles, and binary MSPs with filled circles. Brown bold error bars with 
caps are for MSPs, where dDM^ji, is larger than Cadm- The grey area shows the region where |ADM| < 3 


accounted for by fitting for DM, resulting in an inaccurate DM 
with larger DM uncertainty than without profile variation. The 
importance of taking pulse profile evolution into account was al¬ 
ready demonstrated in Hassall et al. (2012), where they tested 
the accuracy of the dispersion law through wide-hand simulta¬ 
neous observations with LOFAR at 40-190 MHz and the Lovell 
and Effelsberg radio telescopes at 1.4 and 8 GHz, respectively. 

To see if this was indeed the case for our MSPs, we calcu¬ 
lated DM values using both the Tempo2 and pdmp (van Straten 
et al. 2012) programs. For pdmp, DM optimisation is based on 
the signal-to-noise ratio of the entire profile, while for Tempo2 
it is based on cross-correlation between template and observed 
profile. In the ideal case of a simple, unscattered Gaussian pro¬ 
file, identical across the whole observing band, both Tempo2 and 
pdmp should provide the same DM measurements. For scattered 
profiles, one would expect pdmp to give larger DM values due to 
overcompensation for scatter-broadening, although a DM mea¬ 
surement by Tempo2 would be biased as well. If there is a vari¬ 
ation of the profile in amplitude due to changes in the bandpass 
gain or intrinsic to the pulsar itself, then this should not matter 
for pdmp, but for Tempo2 it will depend on how well the match 
is between the template and profiles across the band. Therefore, 


generally in the case of any profile variation one would expect 
to see a discrepancy between DM measurements provided by 
Tempo2 and pdmp. 

Columns 11 and 12 of Table 5 list the Tempo2 and pdmp DM 
measurements for the best individual observations. To charac¬ 
terise the difference between these values we also calculated the 
parameter eoM as; 

_ |I9Mxei;jipQ2 DMpdmpI 

fDM = y , , = > 

V^DM2,^p„ 2 + bDM^d^p 

where bDMTempo 2 and ADMpdmp are DM uncertainties from 
Tempo2 and pdmp. For 30 MSPs in our sample the value of 
fDM is less than one, and for all of them we do not see any 
profile variation. For the other 15 pulsars eoM > F and all of 
them manifest a visible change of profile shape and/or ampli¬ 
tude across the band. Among these 15 pulsars, six MSPs have 
euM > 2, namely PSRs J1640+2224 (eoM = 2.1), J1853-H1303 
(2.2), B1937-H21 (2.65), J2051-0827 (2.65), J1810H-1744 (3.1), 
and J1816-I-4510 (3.4) in ascending value of sdm- Figure 2 al¬ 
ready shows the profile evolution of PSR B1937-1-21, and that 
profile becomes completely scattered out below 150 MHz. Pro- 
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Pulse phase Pulse phase Pulse phase 

Fig. 9. Profile evolution for three MSPs within the LOFAR HBA band. The number on the left in each panel gives the central frequency of the 
corresponding sub-band in MFlz. The profile in grey on top of each panel is the average profile within the 78-MFtz wide band. The profiles for 
each sub-band are normalised by its noise rms deviation. Each sub-band profile is overlaid with the average profile of the total band, in grey. To 
keep the integrated flux densities the same, the overlaid total-band average profile (in grey) is normalised by the ratio of areas under the total-band 
profile and the sub-band profile. 


file evolution for the other three MSPs with the largest eoivi val¬ 
ues is shown on Fig. 9. One can see how the profiles of PSRs 
J1810H-1744 and J1816+4510 are scattered towards lower fre¬ 
quencies. For PSR J2051-0827 scattering does not seem to be 
relevant, but the pulse amplitude differs significantly at the cen¬ 
tre and edges of the band. Therefore, the eDM parameter could be 
used as indirect indication for noticeable profile variation within 
the observing band, and whether it is needed to be taken into 
account for accurate DM measurements. 

It is clear that for accurate measurements of DM, especially 
at low frequencies, profile evolution must be taken into account. 
Even very small profile variations at higher frequencies can 
also systematically bias the derived arrival times and must be 
taken into account for high-precision pulsar timing. Frequency- 
dependent profile templates that take profile broadening due to 
scattering into account must be used for this purpose, as im¬ 
plemented in the software package PulsePortraiture'* (Pennucci 
et al. 2014), or using a Generative pulsar timing analysis frame¬ 
work (Lentati et al. 2015). Recently, Liu et al. (2015) considered 
the double-lens model that reproduces the parameters of the ob¬ 
served diffractive scintillation with high accuracy, and discussed 
its prospects for removing scattering to improve pulsar timing. 

4.5. Multifrequency MSP profiles 

Here we discuss qualitatively the wide-hand (100-1510MHz) 
profile evolution and compare our findings to other available 
low-frequency data. Quantitative comparison including high- 
energy light curves will be presented in subsequent work. 

For qualitative comparison we used radio profiles from the 
EPN Database of Pulsar Profiles*® (Lorimer et al. 1998). Eor 
these EPN profiles timing information is not available, there¬ 
fore we visually aligned them by the phase of the profile peaks. 

https://github.com/permucci/PulsePortraiture 

http://www.epta.eu.org/epndb/ 


In case of two or more components in a profile we followed 
alignment from previously published data (Kramer et al. 1999) 
or our own reasoning. These profiles are presented in Pig. 10, 
where we also plot (where available) the previous low-frequency 
Pushchino profiles at 102/111 MHz from Kuzmin & Losovsky 
(1999). LOPAR profiles (blue) are generally of better quality, 
with higher signal-to-noise ratio and time resolution. More im¬ 
portantly, they do not suffer from intra-channel dispersion smear¬ 
ing, as the Pushchino profiles do. Hence, they provide a clearer 
view of profile evolution at low frequencies. Stappers et al. 
(2008) detected eight MSPs at 115-175 MHz with the WSRT. 
Profiles of two pulsars, PSRs J0034-0534 and J2145-0750, are 
presented in their paper (Stappers et al. 2008), and their shapes 
fully resemble the shape of the LOPAR profiles. Contrary to our 
LOPAR observations in 2013 and BSA observations by Kuzmin 
& Losovsky (1999) in 1999, Stappers et al. (2008) did not detect 
PSR J0218-1-4232 in 2008. This could simply be due to a lack 
of sensitivity of WSRT, but it might also hint in favour of reg¬ 
ular variations in the scattering conditions in the ISM on time 
scales of several years towards this pulsar. Regular monitoring 
with LOPAR will answer this question for this and other pulsars 
in our MSP timing campaign. 

Based on how the MSP profiles in Pig. 10 evolved in the 
LOPAR HBA band, we can clearly divide the presented MSPs 
into two distinct groups: a) with LOPAR profiles strongly af¬ 
fected by scattering, and b) LOPAR profiles that do not show 
profile broadening due to scattering, or do so very weakly. Only 
two MSPs from Pig. 10, PSRs J0621H-1002 and 11911-1114, 
fall into the first group. Prom Pig. 1 we can certainly conclude 
that there are more MSPs with significantly scattered profiles 
in the LOPAR HBA band, e.g. PSRs J0218-H4232, B1937-b21, 
B1953-b29, but we will leave the discussion of their profile evo¬ 
lution for a future paper. Both MSPs with scattered profiles in 
Pig. 10 have relatively large DMs, 37 and 31 pccm“^ for PSRs 
10621 -H1002 and 11911-1114, respectively. They do not seem to 
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Fig. 10. LOFAR MSP average profiles (blue) in comparison with profiles at other frequencies (black) mostly from the EPN Database of Pulsar 
Profiles (Lorimer et al. 1998), where available. Each profile is normalised to peak amplitude. These profiles are taken from: Pushchino - Kuzmin 
& Losovsky (1999) and Kuzmin & Losovsky (1996), WSRT - Stappers (private communication), Effelsberg - Kramer et al. (1998), GB140 - 
Camilo et al. (1996a) and Sayer et al. (1997), Parkes - Bailes et al. (1997) and Lorimer (1994), Arecibo - Camilo et al. (1996a), Foster et al. 
(1993) and Camilo et al. (1993), and Lovell - Stairs et al. (1999), Bailes et al. (1997), Gould & Lyne (1998), and Lorimer et al. (1996). 


show the development of other profile components. In the case 
of PSR J0621 + 1002 the visual alignment was arbitrary as the 
scattering moves the peak of the profile to later pulse phases (see 
e.g. the LBA profile of PSR J0034-0534 in Fig. 5). 

The other 11 MSPs from Fig. 10 form the second group, all 
with unscattered or weakly scattered LOFAR profiles. In fact, for 
PSR J2145-0750 Kuzmin & Losovsky (1999) claimed that the 
separation between components decreased at 102 MHz, which 
they interpreted as an effect caused by the quadrupole magnetic 
field. Kramer et al. (1999) did not see this abnormal frequency 
dependence at higher frequencies neither for this nor the other 
MSPs they studied. This conclusion by Kuzmin & Losovsky 
(1999) might be due to the fact that they did not use coher¬ 
ently dedispersed profiles and they were affected by intrachan¬ 
nel dispersion that made it more difficult to measure the sepa¬ 
ration. Kramer et al. (1999) also argued against this conjecture 


by Kuzmin & Losovsky (1999) because the trailing component 
itself consists of two overlapping components (see also Sallmen 
1998). Similarly to Kuzmin & Losovsky (1999) we do not see 
a decrease in the component separation for PSR J2145-0750 at 
LOFAR frequencies in comparison with higher radio frequen¬ 
cies. Also, there is no broadening of the profile due to scattering 
in the HBA band, and we could only measure the scattering time 
to be less than 0.7 ms in the LBA band at 57.7 MHz. 

It is quite interesting that for those MSPs with overlap¬ 
ping profile components (except maybe for PSR J0034-0534), 
which constitute half of the second group, the leading com¬ 
ponent weakens while the trailing component becomes domi¬ 
nant. In particular, this is clearly evident for PSRs J1022+1001, 
B1257+12, J1640+2224, J2145-0750, and J2317+1439. This 
might explain the apparent narrowing of MSP profiles at low fre¬ 
quencies. Dyks et al. (2010) demonstrated that such an enhance- 
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ment of the trailing component in MSP profiles is caused by 
aberration and retardation. There is also an opposing co-rotation 
effect, namely the weakening of the trailing side of the profile 
due to asymmetry of curvature radiation about the dipole axis. 
Therefore, whether an MSP profile shows such a rotation asym¬ 
metry is determined by the dominant effect. Dyks et al. (2010) 
showed that a sharp edge of the trailing component should be 
accompanied by a maximum in the gradient of the polarisation 
angle curve at the same spin phase. We indeed confirm this for 
PSR J1022-H1001 (see Noutsos et al. 2015). 

The profile evolution of PSR J1024-0719 is a very inter¬ 
esting case. At present we aligned the LOFAR profile with the 
peak of the trailing component of the 436-MHz Parkes profile, 
because i) the tail is reminiscent of the tail in the high-frequency 
profiles, and ii) the apparent tendency for a weakening leading 
component from 1.4 GHz to 436 MHz. However, the Pushchino 
profile at 102 MHz is inconclusive and looks somewhat similar 
to the two-component profile at 436 MHz. Proper profile align¬ 
ment is needed to make a firm conclusion. 

In general, without considering the broadening of profiles 
due to scattering, the MSP profiles continue to show constant 
separation between profile components and their widths at fre¬ 
quencies below 200 MHz. The same lack of profile evolution of 
MSPs was reported by Kramer et al. (1999) between frequen¬ 
cies of several GHz down to about 400 MHz. This is different 
from what is observed for normal pulsars and indicates that the 
emission must originate from a very compact region in the mag¬ 
netosphere. 

5. Summary 

We have carried out a LOFAR census of 75 Galactic MSPs in 
the frequency range 110-188 MHz, and detected 48 of them (25 
of which were detected at these frequencies for the first time). 
We have also detected three MSPs out of nine observed with 
the LOFAR LBAs in the frequency range 38-77 MHz, namely 
PSRs 10030-1-0451, J0034-0534, and J2145-0750. For the de¬ 
tected MSPs; 

- We provide average profiles for all detected MSPs, with 
about half of them being the best high-quality profiles so 
far at these frequencies. For 25 MSPs the presented pro¬ 
files are the first at 110-188 MHz. About 35% of the MSPs 
show strong narrow profiles, another 25% exhibit undoubt¬ 
edly scattered profiles, and the rest have low signal-to-noise 
profiles. 

- We measure the mean flux density at 110-188 MHz and 
compare it with the predicted values derived from high- 
frequency observations. For at least a third of the MSPs, 
the main uncertainty on the predicted flux density originates 
from poor knowledge of the spectral index. For about an¬ 
other third, the measured flux densities are still lower than 
the predicted ones even within the uncertainty on their spec¬ 
tral indices. In Sect. 3.2.1 and Sect. 3.2.3 we consider differ¬ 
ent factors that could affect our flux measurements. 

- We also measure the effective pulse width at 150 MHz. For 
the majority of the detected MSPs, their pulse duty cycle, 
6, is less than 20%. Almost all MSPs with 6 > 20% show 
profiles widened by scattering. 

- There is no clear dependence of the MSP detectability on 
the predicted scattering. We discuss this and other possible 
factors that affect the MSP detectability in Sect. 4.1. 

- We present average values of DM and their offsets from the 
catalogue values. For 14 MSPs in our sample the absolute 


value of their DM offsets is more than three times larger than 
their errors, and we do not see any trend of DM offsets be¬ 
coming larger for larger epoch offsets. 

- Finally, we qualitatively compare LOFAR profiles for a sub¬ 
set of MSPs with their profiles at higher radio frequencies. 
If not broadened by scattering, they show apparent consis¬ 
tency in the width of their profile components and their sep¬ 
aration within the frequency range from 150 MHz to a few 
GHz. This is very different from what is observed for normal 
pulsars and suggests a compact emission region in the MSP 
magnetosphere. 
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